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ABSTRACT 


A  STUDY  OF  A  METEOR  BURST  COMMUNICATION  SYSTEM. 


This  study  discusses  the  use  of  the  meteor  burst  channel 
in  communication  systems  and  has  three  main  thursts: 

1)  The  development  and  enhancement  of  an  accurate  and 
reliable  channel  model  based  on  recently  available  empririal 
data.  Analysis  of  this  model  results  in  analytical 
expressions  for  communication  parameters  such  as  channel 
duration  and  throughput  to  be  used  as  design  and  analysis 
tools. 

2)  The  optimization  of  throughput  for  fixed  transmis¬ 
sion  rate  under  the  constraint  of  a  given  maximum  bit  error 
rate.  This  will  demonstrate  the  room  for  improvement  in 
existing  systems  using  constant  transmission  rate. 

3)  The  feasibility  of  efficiently  communicating  over 
the  MBC  using  variable  bit  rate  and  employing  a  feedback 
protocol  to  monitor  the  channel.  This  approach  will 
dramatically  improve  the  throughput  in  comparison  with 
constant  transmission  rate  systems. 
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1  STATEMENT  OF  THE  PROBLEM  AND  REVIEW 


There  is  currently  an  urgent  need  for  alternate 
channels  for  communication.  The  large  variety  and 
quantity  of  users  throughout  the  world  have  caused  all  of 
the  existing  channels  to  become  extremely  congested.  A 
world  wide  search  for  such  alternate  channels  is  now 
being  conducted  for  both  present  day  and  future 
communication  systems  through  the  use  of  different  types 

of  channels  such  as  optical,  laser,  cable,  satellite, 

( 

etc ^14-/  IS,"  17] 

f 

-  The  particular  application  of  Beyond-Line-of-Sight 
(BIOS)  communications,  which  uses  the  High  Frequency  (HF) 
spectral  band  from  3  to  30  MHz,  is  very  important  to  many 
current  and  potential  users.  However,  HF  is  very 
sensitive  to  solar  disturbances  such  as  sunspot  activity, 
solar  storms  and  other  galactic  phenomena,  as  well  as 
multipath  returns  from  both  ground  and  atmosphere, 

weather  related  attenuation,  and  other  degrading  factors 

c.~ 

[18,19],  Because  of  these  factors  a  channel  for 
communication  is  needed  which  has  low  congestion,  is 


robust  and  relatively  indestructible,  has  very  little 
outside  interference  and  does  not  compete  for  bandwidth 
with  existing  communication  systems.  A>_The  research 
proposed  herein  represents  such  an  alternative  mode  for 
communication.  _ _ 

Meteors  are  small  aggregates  of  matter  which  upon 
entering  the  Earth's  atmosphere  burn  up  and  form  long 
columns  or  trails  of  ionized  particles  in  the  upper 
atmosphere  at  altitudes  of  80  to  120  km.  The  ionized 
columns  diffuse  within  seconds  yet  are  able  to  support 
radio  communication  over  a  range  of  distance  up  to  2000 
km  by  reflecting  a  radio  signal  and  thus  provide  us  with 
the  so-called  meteor  burst  channel  (MBC) .  The  frequency 
range  used  is  from  30  to  100  MHz.  The  lower  limit  is 
dictated  by  the  need  to  be  above  the  HF  range  and  avoid 
ionospherically  reflected  signals.  The  upper  limit  is  set 
primarily  by  equipment  sensitivity  limitations  since 
reflections  at  the  higher  frequencies  are  weaker  than 
those  at  lower  frequencies. 

Excluding  showers  (annual  Pleiades  shower)  the 
occurrence  of  a  meteor  burst  appears  in  time  to  be  random 
with  a  Poisson  arrival  and  an  average  interval  time  on 


the  order  of  several  seconds  to  minutes  depending  upon 
time  of  day,  season,  and  global  location,  and 
communication  system  parameters. 

Two  types  of  trails  are  described  in  the  literature: 
Underdense  trails  which  are  more  frequent  and  have  an 
electron  line  density  below  1014  e/m;  and  overdense 
trails  which  are  less  frequent  and  have  electron  line 
densities  of  more  than  1014  e/m.  The  power  received  from 
an  underdense  trail  is  approximated  as  an  exponential 
time  function  whose  initial  amplitude  is  proportional  to 
the  square  of  the  electron  line  density  and  the 
exponential  decay  time  constant  varies  randomly  from  one 
trail  to  the  next.  The  power  received  from  an  overdense 
trail  is  modeled  by  a  complex  time  function  whose  shape 
for  a  given  set  of  communication  system  parameters  is 
determined  by  the  value  of  the  electron  line  density. 

In  the  late  fifties  and  early  sixties,  a  considerable 
amount  of  research  was  conducted  on  the  use  of  meteor 
burst  communications  for  beyond  line-of-sight  (BLOS) 
transmission  of  digital  data.  In  recent  years,  interest 
in  meteor  burst  has  been  renewed  because  of  several 
factors.  The  development  of  microprocessors  (enabling 
inexpensive  implementation  of  sophisticated  system 
control  procedures)  and  inexpensive  solid-state  memories 
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(for  the  buffering  required  to  interface  a  burst 
transmission  scheme  to  constant  data  sources  and  sinks) 
is  partly  responsible.  In  addition,  the  nuclear 
survivability  of  the  meteor  burst  medium  is  superior  to 
other  BLOS  media  such  as  satellite  and  HF  skywave, 
resulting  in  considerable  interest  in  the  military 
community . 


1.2  Identification  of  the  Problem 


The  research  proposed  herein  represents  such  an 
alternative  mode  for  BLOS  communication.  It  relates  a 
novel  approach  for  BLOS  communication.  The  approach  is 
aimed  at  efficaciously  using  a  previously  under  utilized 
and  inefficiently  used  method  of  communication,  namely, 
the  Meteor  Burst  Channel  (MBC)  [1,  9].  This  medium 

possesses  all  of  the  qualities  mentioned  above  for  a  new 
communications  channel.  The  channel  operates  in  the 
relatively  unused  lower  portion  of  the  Very  High 
Frequency  (VHF)  band  ranging  from  30  to  100  MHz.  It 
avoids  the  degradations  exhibited  by  HF  on  the  low  end 
and  ionospheric  phenomena  at  the  high  end.  This  region 
of  the  spectrum  provides  a  sufficiently  large  bandwidth 
for  efficient  data  communication.  It  is  not  easily 
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destroyed.  It  has  a  privacy  feature  inherent  in  its 
structure  [1,9]  and  by  virtue  of  its  NnevnessN  does  not 
interfere  with  existing  systems.  For  MBC,  the  antennas 
have  a  high  elevation  angle  and  rarely  interfere  with 
present  use  of  the  frequency  band  when  such  is  present. 

This  study  seeks  to  define  systematically  and 
analytically  the  limits  of  transmission  of  digital 
information  over  the  meteor  burst  channel  by  using  an 
improved  model  to  develop  the  mathematical  relationships 
among  link  geometry,  communication  circuitry  parameters, 
physical  meteor  characteristics,  modulation  technique  and 
a  given  communication  protocol.  Furthermore,  the  re¬ 
search  will  be  aimed  at  the  determination  of  optimal 
average  throughput  for  both  constant  and  variable 
transmission  rates  under  the  employment  of  an  information 
feedback  scheme  for  reliable  channel  monitoring. 

These  results  will  serve  as  a  vital  design  and 
analysis  tool  for  the  communication  system  engineer  using 
the  MBC.  Furthermore,  they  will  demonstrate  the  viability 
and  potential  improvement  over  past  and  currently 
implemented  systems.  Finally,  they  will  identify  the 
necessary  and  challenging  areas  for  future  research. 


Meteor  Burst  communications  has  vide  ranging  effects 
for  modem  society.  It  offers  simple  and  reliable 
communications  and  can  operate  in  the  simplex,  half-du¬ 
plex  and  full-duplex  mode  over  distances  of  200  to  1200 
miles.  [9,22].  Efficient  use  of  this  channel  will 
result  in  increased  data  throughput  for  applications  such 
as  facsimile  transmission  (FAX) ,  TELEX,  transmission  of 
computer  data,  radio  amateur  usage,  voice  burst 
transmission  using  phonomes  (i.e.,  transmission  in  which 
the  voice  is  digitally  encoded  using  rates  of  50  to  250 
bits  per  second)  and  military  communication  systems. 
These  types  of  communications  and  data  transmission 
facilities  involve  transactions  of  banks,  universities 
and  medical  centers,  corporate  business  concerns  such  as 
online  systems,  and  government  agencies,  in  addition, 
recreational  and  pleasure  or  leisure  activities  could 
make  use  of  this  mode  of  communication. 

If  shown  to  be  feasible  and  commercially  cost 
effective  business  as  well  as  home  users  could  provide  a 
potentially  large  market.  The  advent  of  the  microproces¬ 
sor  will  reduce  the  cost  of  such  systems  thus  enabling 
inexpensive  implementation  of  the  control  requirements  to 
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be  met  [9,22,23].  This  potential  for  low  cost  will  open 
up  the  home  and  business  market  which  has  already 
invested  millions  of  dollars  in  personal  computers  and 
modems.  Furthermore,  the  number  of  these  types  of  users 
is  increasing  daily. 

Some  of  the  government  agencies  interested  in 
alternate  communications  means  are  the  Department  of 
Agriculture,  the  Department  of  Energy  (DOE),  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA) ,  the 
Department  of  Defense  (DOD) ,  the  National  Aeronautics  and 
Space  Administration  (NASA)  and  the  Defense  Communica¬ 
tions  Agency- ( DC A)  [9,22]. 


rsica] 


rties  of  Meteor  Trails 


Nearly  all  of  the  present  knowledge  regarding  the 
physical  properties  of  meteors  has  been  obtained  from 
visual,  optical,  and  radio  observations  of  the  trails 
formed  by  the  meteoric  particles  as  they  enter  the 
atmosphere  of  the  earth.  In  this  section  the  properties 
of  the  meteoric  particles  themselves  and  the  trails  which 
they  form  will  be  reviewed.  The  mathematical  model  and 
conceptual  framework  reflecting  physical  and  radio 
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properties  of  the  trail  will  be  presented  in  chapter  2. 
The  review  of  meteors  in  sections  1.4.1-  1.4.4  is  taken 
directly  from  Sugar's  [1]  1964  review  paper. 


1.4.1  Meteoric  Particles 


For  the  purposes  of  this  discussion  the  term 
"Meteors"  will  be  used  to  apply  to  those  particles 
entering  the  earth ' s  atmosphere  that  are  burned  by 
frictional  heating.  This  definition  includes  the  very 
small  particles,  the  micrometeorites,  that  slowly 
settle  through  the  atmosphere  without  being  destroyed. 
It  also  includes  the  large  meteors  which  manifest 
themselves  as  fireballs  or  the  even  larger  ones  which 
reach  the  earth's  surface  as  meteorites.  The  microme¬ 
teorites  are  not  of  concern  because  even  though  they 
are  the  most  numerous  of  the  various  types,  they  enter 
the  atmosphere  too  slowly  to  produce  any  significant 
ionization.  The  large  meteoroids,  although  they 
produce  substantial  ionization,  are  of  little  concern 
here  because  their  rate  of  occurrence  is  extremely  low. 

Some  of  the  properties  of  meteoric  particles  are 
summarized  in  Table  1  [24],  [25].  The  particles  of 
interest  in  meteor  propagation  are  those  with  masses  in 


tha  rang*  103  to  10“7  g  and  dimensions  in  th*  range  8 
cm  to  40  microns.  Before  being  trapped  by  the 
gravitational  field  of  the  earth  theae  particlea  move 
in  orbita  around  the  aun.  Their  conpoaition  ia 
uncertain;  however/  they  appear  to  be  alnoat  entirely 
of  coaetary  origin.  A  aubatantial  fraction  of  them  are 
not  aingle  solid  particles  but  fragile  loosely-bound 
agglomerates/  sometimes  called  •’dustballB", 

The  meteors  can  be  divided  into  two  classes/  the 
shower  meteors  and  the  sporadic  meteors.  The  shower 
meteors  are  collections  of  particles  all  moving  at  the 
same  velocity  in  fairly  well-defined  orbits  or  streams 
around  the  sun.  Their  orbits  intersect  the  orbit  of 
the  earth  at  a  specific  time  each  year  and  at  these 
times  the  well-known  meteor  showers  are  observed.  In 
the  cases  where  the  particles  are  uniformly  distributed 
around  the  orbit,  the  size  of  the  shower  varies  little 
from  year  to  year.  If,  on  the  other  hand,  there  is  a 
concentration  of  particles  within  the  orbit,  the  extent 
of  a  shower  can  vary  substantially  in  successive  years. 
(The  Leonid  shower  of  1833  is  an  example  of  an  unusual 
display  of  meteor  activity.  During  the  peak  of  the 
shower  as  many  as  20  meteor  trails  were  often  seen  by  a 
single  observer  in  one  second  .)  The  shower  meteors, 
while  the  most  spectacular,  account  for  only  a  small 


fraction  of  all  meteors.  It  is  the  nonshower  or 
sporadic  meteors  that  comprise  nearly  all  of  the 


meteors 

of 

interest  in 

radio 

propagation.  These 

meteors 

are 

those 

which 

do  not 

have  well-defined 

streams 

but 

rather 

seem 

to  move 

in  random  orbits . 

Thus,  whereas  shower  meteors  appear  to  be  coming  from  a 
specific  point  in  the  sky — the  radiant  point  for  the 
shower — sporadic  meteors  have  radiants  that  appear  to 
be  randomly  distributed  over  the  sky. 

The  relation  between  sporadic  and  shower  meteors  is 
not  clear  at  the  present  time.  It  is  possible  that  the 
sporadic  meteors  represent  the  final  stages  in  the 
decay  of  meteor  streams.  On  the  other  hand  there  is 
some  indication  that  the  sporadic  meteors  are  in  fact 
distributed  in  a  large  number  of  relatively  small 
groups  and  that  the  earth  is  immersed  in  about  11  such 
groups  at  a  time.  Further  work  is  needed  to  resolve 
the  origins  of  the  sporadic  meteors  [1]. 

Shower  meteors  are  most  easily  recognized  in  terms 
of  their  radiants,  velocities,  and  time  of  occurrence, 
since  these  parameters  are  relatively  fixed.  The 
radiants  and  times  of  occurrence  of  sporadic  meteors 
are  random.  Their  radiant  points  are  not,  however, 
uniformly  distributed  in  the  sky  but  are  for  the  most 


part  concentrated  toward  the  ecliptic  plane  (the  plane 
of  the  earth's  orbit)  and  move  in  the  same  direction 
around  the  sun  as  the  earth  moves  [26].  The  orbits  are 
not  uniformly  distributed  along  the  earth's  orbit  but 
are  concentrated  so  as  to  produce  a  maximum  incidence 
of  meteors  at  the  earth  in  July  and  a  minimum  in 
February.  This  variation  in  the  space  density  of 
meteors  is  shown  in  Fig.  1  which  has  been  adapted  from 
the  data  of  Hawkins  [26]. 

The  rate  of  incidence  of  sporadic  meteors  at  the 
earth  is  further  modified  by  two  factors.  The  first  of 
these — resulting  in  a  regular  diurnal  variation  in 
meteor  rate  is  illustrated  in  Fig.  2  .  On  the  morning 
side  of  the  earth,  meteors  are  swept  up  by  the  forward 
motion  of  the  earth  in  its  motion  around  the  sun.  On 
the  evening  side  the  only  meteors  reaching  the  earth 
are  those  which  overtake  it.  This  results  in  a  maximum 
occurrence  rate  around  6  a.m.  and  a  minimum  rate  around 
6  p.m.  The  ratio  of  maximum  to  minimum  depends  on  the 
latitude  of  the  observer.  A  further  minor  seasonal 
variation  is  introduced  because  of  the  tilt  of  the 
earth's  axis  relative  to  the  ecliptic  plane  [27].  This 
variation,  also  dependent  on  the  latitude  of  the 
observer,  can  change  the  expected  hourly  rates  by 
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The  Base  distribution  of  sporadic  meteors  is  such 
that  there  are  approximately  equal  total  masses  of  each 
size  of  particle.  There  are,  for  example,  10  times  as 
many  particles  of  mass  10~4  g  as  there  are  particles 
of  mass  10“3  g.  This  approximate  relation  between 
particle  mass  and  number  is  given  in  Table  1.  The  mass 
distribution  of  shower  meteors  is  somewhat  similar  to 
that  for  the  sporadic  meteors  with  the  important 
difference  that  there  are  more  large  particles  relative 
to  the  number  of  smaller  ones  than  for  the  sporadic 
meteors  . 

The  velocities  of  meteors  approaching  the  earth  are 
in  the  range  11.3  to  72  km/sec.  The  lower  limit  is  the 
escape  velocity  for  a  particle  leaving  the  earth  and  is 
therefore  the  lowest  velocity  that  a  particle  falling 
toward  the  earth  can  have.  The  upper  limit  is  the  sum 
of  two  components,  a  30  km/sec  component  associated 
with  the  earth  in  its  orbit  around  the  sun,  and  a  42 
km/sec  component  associated  with  the  meteor  itself. 
This  latter  velocity  is  the  escape  velocity  for  a 
particle  leaving  the  solar  system.  Nearly  all 
observations  have  indicated  that  meteor  velocities  fall 
in  the  above  range  and  that  the  meteors  are  thus 
members  of  the  solar  system. 
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1)  Formation:  As  a  meteoric  particle  enters  the 
earth's  atmosphere  it  collides  with  air  molecules  which 
then  become  trapped  in  its  surface.  The  impact  energy 
produces  heat  which  evaporates  atoms  from  the  meteor 
and  these  move  off  with  a  velocity  which  is 
substantially  equal  to  that  of  the  meteor.  Collisions 
between  these  high  velocity  atoms  and  the  surrounding 
air  result  in  the  production  of  heat,  light,  and 
ionization,  distributed  in  the  form  of  a  long,  thin, 
paraboloid  of  revolution  with  the  particle  at  the  head 
of  it.  The  electron  line  density  in  the  trail  is 
proportional  to  the  mass  of  the  particle.  In  the 
evaporation  process  each  original  impact  frees  many 
meteoric  atoms  and  thus  the  total  mass  of  air  molecules 
striking  the  meteor  is  small  compared  to  the  meteoric 
mass  itself.  As  a  consequence  the  velocity  of  the 
meteor  remains  quite  constant  until  the  meteor  is 
nearly  completely  evaporated.  (Those  meteors  that 
appear  to  be  dustballs  rather  than  solid  particles 
break  up  or  "fragment"  upon  entering  the  atmosphere  and 
proceed  as  a  group  of  smaller  particles.  For  this  case 


the  above  description  appears  to  be  correct  if  it  is 
applied  to  the  individual  particles  after  fragmentation 
ocurs . ) 

2)  Heights:  As  a  meteoric  particle  approaches  the 
earth  no  appreciable  ionization  if  formed  until  the 
particle  enters  the  relatively  dense  air  at  heights 
below  about  120  km.  Above  that  height  collisions  of 
the  particles  with  air  molecules  are  not  frequent 
enough  to  be  of  significance.  As  the  particle 
traverses  the  region  below  about  120  km  it  vaporizes 
rapidly  and  most  particles  are  completely  evaporated 
before  reaching  80  km.  The  relatively  small  thickness 
of  the  meteor  region  is  a  result  of  the  rapid  change  in 
air  density  which  occurs  over  the  height  range  quoted. 
At  120  km  the  mean  free  path  is  5.4  m  and  this  decreaes 
in  an  approximately  logarithmic  manner  so  that  at  80  km 
it  is  only  3.8  mm. 

The  height  distribution  of  trails  varies  somewhat 
with  particle  characteristics.  The  higher  velocity 
particles  produce  trails  at  the  higher  heights  with  the 
mean  trail  height  increasing  to  about  10  km  as  velocity 
increases  from  less  than  15  km/sec  to  greater  than  60 
km/sec  [1].  The  particles  of  higher  mass  produce 
maximum  trail  ionization  at  lower  heights.  Over  the 


Bass  range  of  103  to  10~7  g  the  height  variation  is 
around  44  km.  There  is  also  a  height  variation  with 
the  zenith  angle  of  the  trail  orientation,  with  the 
larger  angles  corresponding  to  greater  heights.  A 
variation  of  about  13  km  is  associated  with  nearly  the 
whole  range  of  zenith  angles. 

3)  Lengths:  The  lengths  of  meteor  trails  are 
primarily  dependent  on  particle  mass  and  zenith  angle. 
Typical  lengths  range  up  to  50  km  with  the  most 
probable  trail  length  for  sporadic  meteors  being  15  km. 
(Several  definitions  of  "length"  can  be  used  and  the 
one  chosen  here  uses  as  end  points  the  points  with  a 
threshold  value  chosen.) 

4)  Initial  Radius:  Until  recently  it  has  been 
assumed  that,  at  the  time  of  their  formation,  trails 
had  an  initial  radius  of  the  order  of  the  mean  free 
path  or  at  most  about  14  times  this  radius.  However, 
photographic  and  radio  measurements  have  suggested  that 
the  initial  radius  is  significantly  larger  than 
indicated  above  and  that  this  increase  is  probably 
associated  with  the  dustball  and  fragmentation 
hypothesis. 
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The  trail  radii  indicated  by  these  measurements  are 
in  the  range  o  to  1.2  m  (with  a  mean  value  of  0.65  m) 
for  the  photographic  work,  and  0.55  to  4.35  m  for  the 
radio  work.  The  0.55  m  was  for  a  height  of  81  km  where 


the  mean  free  path  is  only  5  x  10  m.  The  results  at 

121  km  are  more  nearly  consistent  with  Manning's 

hypothesis  since  the  initial  radius  observed  is  4.35  m 
and  the  mean  free  path  is  5  m. 

5)  Dissipation:  once  the  trail  is  formed  it 

expands  by  diffusion  [21]  at  a  relatively  low  rate, 

producing  a  radial  distribution  of  material  that  is 
approximately  Gaussian.  The  quantity  (4  Dt  +  r  )  may 
be  taken  as  the  approximate  radius  of  the  trail  after  a 
time  t  where  D  is  the  diffusion  coefficient  and  r  is 
the  initial  radius  of  the  trail.  D  varies  from  1 
m2/sec  at  an  85  km  height  to  140  m2/sec  at  115  km. 
Thus  after  one  second  a  trail  will  have  a  radius  in  the 
range  2  to  20  m. 

The  practical  lifetime  of  a  trail  is  of  course 
dependent  on  the  means  for  detecting  it.  Most  trails 
detected  by  radio  means  are  those  resulting  from  small 
particles  and  these  last  only  a  fraction  of  a  second. 
The  larger  particles  produce  more  densely  ionized 
trails  and  trails  with  durations  of  the  order  of  a 


minute  are  observed  several  times  per  day.  Trails 
with  durations  of  the  order  of  an  hour  or  more  are 
extremely  rare. 


The  dissipation  of  trails  is  further  complicated  by 
the  presence  of  winds  in  the  meteor  region.  At  the 
time  of  formation  trails  are  quite  straight  but  they 
are  rapidly  deformed  by  these  winds  which  have  typical 
velocities  of  the  order  of  25  m/sec  and  vertical 
gradients  with  mean  relative  maxima  of  about  100 
m/sec/km.  A  wind  shear  of  this  magnitude  can  rotate 
part  of  a  trail  through  an  angle  of  5°  in  one  second. 


1*4.3  Reflection  Properties  of  Trails 

The  distribution  of  energy  reflected  by  a  meteor 
trail  is  a  function  of  many  variables.  The  ionization 
density  distribution  across  and  along  the  trail,  the 
orientation  of  the  trail,  the  radio  wavelength,  the 
polarization  of  the  incident  wave  relative  to  the 
trail,  motion  of  the  trail  either  as  part  of  the 
process  of  formation  or  due  to  ionospheric  winds,  and 
the  straightness  of  the  trail  are  all  significant.  In 
discussing  the  reflection  properties  it  is  convenient 
to  divide  the  trails  into  two  classes  underdense  trails 


and  overdense  trails  and  to  examine  the  properties  of 
each  class  independently.  Underdense  trails  are  those 
wherein  the  electron  density  is  low  enough  so  that  the 
incident  wave  passes  through  the  trail  and  the  trail 
can  be  considered  as  an  array  of  independent 
scutterers.  Overdense  trails  are  those  wherein  the 
electron  density  is  high  enough  to  prevent  complete 
penetration  of  the  incident  wave  and  to  cause 
reflection  of  waves  in  the  same  sense  that  the  ordinary 
ionospheric  reflections  occurs.  A  rough  sorting  of 
trails  into  these  two  categories  can  be  done  on  the 
basis  of  trail  lifetime  or  duration.  At  long 
wavelengths  the  underdense  trails  have  durations  of 
less  than  about  one  second  while  the  overdense  trails 
have  longer  durations.  At  long  wavelengths  the 
effective  duration  of  a  trail  is  large  compared  to  the 
time  it  takes  the  trail  to  form,  and  the  trail  may  be 
considered  to  have  a  cylindrical  shape. 

The  solutions  to  be  considered  here  are  useful 
approximations  to  the  physical  problems.  They  will 
apply  quite  well  to  some  of  the  trails  observed  and 
rather  poorly  to  others.  A  complete  analysis  of  the 
reflections  from  even  a  relatively  simple  trail  would 
be  far  too  complex  to  be  of  any  practical  use  or 
interest. 
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1.4.4 


In  this  ssction  sobs  additional  practical  aspects 
will  be  discussed  fron  a  qualitative  viewpoint.  The 
discussion  is  directed  toward  the  long  wavelength  cases 
since  relatively  few  short  wavelength  observations  have 
been  made  and  analyzed. 


1)  Long  Wavelength  Reflections  during  Trail 
Formation:  The  transient  state  associated  with  trail 
formation  is  of  interest  since  it  accounts  for  some  of 
the  observed  characteristics  of  trail  reflections.  As 
a  trail  is  being  formed,  but  before  the  meteor  reaches 
the  first  Fresnel  zone,  a  weak  reflection  is  obtained 
from  the  incomplete  trail.  This  comes  primarily  from 
the  part  of  the  trail  corresponding  to  the  shortest 
transmission  path  at  the  instant;  and  when,  as  is 
usually  the  case,  this  is  the  head  of  the  trail,  the 
reflected  signals  are  shifted  in  frequency  because  of 
the  motion  of  the  effective  reflecting  point.  As  the 
meteor  approaches  the  first  Fresnel  zone  for  the  trail 
this  frequency  shift  approaches  zero,  and  thus  the 
received  frequency  decreases  with  time.  The  observed 
frequency  will  of  course  depend  on  trail  orientation, 
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meteor  velocity  and  observed  wavelength.  A  maximum 
shift  of  the  order  of  5  kHz  is  possible  for  50  MHz 
forward-scatter  observations  over  a  1000  km  path. 

2)  Trail  Drift  and  Distortion:  The  effects  of 
ionospheric  winds  are  appreciable  for  trails  which  last 
for  the  order  of  a  second  or  more.  A  Doppler  shift  of 
the  received  frequency  will  be  associated  with  the 
average  wind  velocity  of  the  trail.  For  a  velocity  of 
25  m/sec  this  "body  Doppler"  can  be  as  large  as  18  Hz 
for  backscatter  observations  at  50  MHz  and  will  be 
somewhat  less  for  forward-scatter  observations. 

The  trail  distortion  resulting  from  wind  shears  can 
lead  to  the  formation  of  several  local  first  Fresnel 
zones  for  the  trail  since  a  distorted  trail  can  have 
several  points  where  the  transmission  path  length  has  a 
local  minimum.  These  local  minima,  or  "glints"  as  they 
have  been  called  [41],  are  strong  scatters  and  the 
received  signal  is  a  composite  of  their  contributions. 
Since  they  are  moving  at  different  velocities  the 
signals  from  each  have  different  Doppler  shifts  and  the 
resultant  composite  signal  fades  in  an  irregular  manner 
[24].  (At  50  MHz  the  fading  rate  observed  for  forward 
scatter  are  of  the  order  of  l-io  Hz.)  In  addition  to 
producing  the  fading  observed  for  long-enduring  trails 


the  wind  shears  can  rotate  a  trail  sufficiently  to 
produce  reflections  when  the  initial  orientation  has 
not  been  suitable.  Thus  these  trails  lose  their  aspect 
sensitivity  as  tine  goes  on,  and  if  their  life  is  of 
the  order  of  10  seconds  they  will  scatter  in  all 
directions . 

3)  Polarization  Effects:  Thus  far  the  only 
indication  that  the  polarization  of  the  incident  wave 
can  be  of  significance  was  the  inclusion  of  the  sin2  (a) 
term  in  each  of  the  transnission  equations  for  the 
forward-scatter  cases.  This  term  accounts  for  the 
foreshortening  of  the  electric  vector  of  the  incident 
wave  which  occurs  when  it  is  viewed  from  the  receiver. 
In  addition  to  this,  a  kind  of  electron  reasonance  can 
occur  which  increases  the  reflection  coefficient  for 
the  trail.  This  is  associated  with  the  restoring  force 
that  the  electrons  experience  as  the  incident  electric 
field  displaces  them  from  their  equilibrium  positions 
in  a  direction  normal  to  the  trail.  This  reasonance 
can  occur  only  for  underdense  trails  and  only  then  for 
trail  diameters  much  less  than  the  wavelength.  This 
reasonance  at  most  doubles  the  reflection  coefficient 
and  thus  can  increase  the  received  power  by  a  factor  of 
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4)  Diversity  Effects:  Underdense  trails  act  as 
small  coherent  sources ,  and  therefore  CW  signals 
scattered  by  them  exhibit  good  space  and  frequency 
correlation.  Correspondingly,  very  little  pulse  broad¬ 
casting  is  observed  when  pulse  signals  are  used.  The 
limited  available  data,  suggest  that  for  a  for¬ 

ward-scatter  path  about  1000  km  in  length,  operating 
near  50  KHz  the  correlation  coefficent  observed  for 
single  underdense  trails  will  fall  to  0.5  for  antenna 
spacings  of  the  order  of  150  km  along  the  path,  for 
antenna  spacings  of  the  order  of  30  km  across  the  path, 
and  for  frequency  separations  greater  than  5  MHz. 

Three  micro-sec  pulses  show  no  appreciable  broadening 
under  these  conditions. 

Overdense  trails,  in  contrast  to  the  underdense 
trails,  tend  to  act  like  relatively  large  sources  and 
therefore  exhibit  much  poorer  space  and  frequency 
correlation  properties.  Again,  as  in  the  underdense 
case,  specific  data  on  these  properties  are  not 

available.  Present  results  suggest  that  the  correla¬ 
tion  observed  for  single  overdense  trails  will  fall  to 
0.5  for  antenna  spacings  of  the  order  of  50  [41]. 

Measurements  using  short  pulses  over  a  1000-km  path 

indicate  that  the  received  composite  signal  can  have  a 
total  time  spread  as  large  as  10  micro-  sec  [44]. 


There  are  three  potential  modes  of  operation  for  a 
meteor  burst  system:  point-to-point,  netted,  and 
broadcast.  While  a  meteor  burst  system  could  be  operated 
in  any  or  all  of  these  modes,  all  known  meteor  burst 
systems  that  have  currently  been  implemented  have  been 
designed  for  point-to-point  applications.  The  general 
discussion  in  the  following  sections  1.5. 1-1. 5. 3  is  taken 
directly  from  Oetting's  [9]  1980  review  paper. 

1.5.1  Point  to  Point  Mode 


The  point-to-point  mode  is,  of  course,  straightfor¬ 
wardly  implemented.  The  only  requisite  for  effective 
system  control  is  the  ability  of  the  transmitting 
terminal  to  discern,  as  accurately  as  possible,  the 
beginning  and  end  of  a  useful  burst.  If  transmission 
begins  too  late  or  ends  too  soon,  valuable  burst  time 
will  be  wasted.  If  the  transmission  extends  past  the 
useful  portion  of  the  burst,  a  high  error  rate  will 
result. 
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In  all  point-to-point  systems,  we  assume  that  a 
feedback  path  Is  available.  In  the  half-duplex  case, 
the  forward  and  return  links  share  the  same  frequency, 
while  in  the  full-duplex  case,  they  use  two  different 
frequencies.  For  achieving  data  transmission  between 
point  A  and  point  B,  perhaps  the  simplest  strategy  is 
to  assign  a  different  frequency  to  each  transmitter  and 
to  continuously  transmit  a  probe  signal  from  the  master 
station  (say,  point  A) .  When  point  B  detects  the 
probe,  he  transmits  a  preamble  followed  by  his 
information.  Point  A  uses  the  preamble  to  synchronize 
his  receiver  and  is  subsequently  able  to  receive  data. 
It  is  usually  required  that  A  acknowledge  the  reception 
of  the  message  or  block  of  data. 

Several  variations  of  this  basic  point-to-point 
approach  have  been  developed.  Early  meteor  burst 
systems  used  a  signal  amplitude  threshold  to  determine 
if  the  probe  signal  were  being  received.  The  JANET 
system  [1]-[13J  improved  upon  this  procedure  by 
substituting  a  signal-to-noise  ratio  measurement  for 
the  fixed  threshold  test,  thereby  reducing  the  false 
alarm  rate.  After  extensive  testing  of  the  JANET  B 
system,  the  SHAPE  Technical  Centre  concluded  that,  even 
with  the  SNR  approach,  a  suitable  compromise  between 
efficient  use  of  the  channel  and  low  error  rate  was  not 


achievable.  Thus,  they  developed  their  own  approach 
using  ARQ  with  one  teletype  character  per  seven-bit 
block.  This  apporach  is  referred  to  as  the  COMET 
system  [5]  and  [9]. 

\ 

The  possibilities  for  netted  meteor  burst  systems 
have  not  been  well  explored  to  date.  The  only  known 
system  capable  of  netted  operation  is  the  Western  Union 
hardware  currently  being  installed  for  the  Department 
of  Agriculture  [28].  This  system  will  consist  of  511 
remote  sites  that  communicate  with  two  master  stations. 
Although  the  Department  of  Agriculture  application 
involves  one-way  transmission  of  data  from  the  remote 
to  the  master  stations,  the  system  can  be  configured  to 
provide  two-way  data  transmission.  A  network  could 
then  be  achieved  by  having  a  master  station  operate  in 
a  s t or e-and- forward  mode. 


1.5.2  LOW  probability  of  Intercept  (LPI)  and 
Antiiam  (AJ)  Consideration 


One  of  the  primary  attractions  of  meteor  burst 
communications  is  its  inherent  privacy  resulting  from 
the  restricted  footprint  of  reflections  from  meteor 
trails.  This  property  has  received  considerable 
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attention  in  tha  meteor  burst  literature  and  has 
prompted  a  large  number  of  industry  proposals  to 
provide  LPI  meteor  burst  systems  for  numerous 
Department  of  Defense  (DoD)  and  non-DoD  agencies. 

Signals  propagating  between  a  specific  transmitter 
and  receiver  can  normally  be  detected  by  an 
eavesdropping  receiver  only  if  the  latter  is  located  in 
an  elliptically  bounded  "footprint"  in  the  vicinity  of 
the  primary  receiver.  By  reciprocity,  interfering 
signals  from  a  remote  transmitter  will  propagate  to  a 
meteor  burst  receiver  only  if  the  interfering 
transmitter  is  within  the  footprint  of  the  desired 
transmitter.  Measurements  indicate  that  the  major  and 
minor  axes  of  the  footprint  are  on  the  order  of  2000 
and  25  mi,  respectively.  Clearly,  the  footprint  is 
not  a  well-defined  boundary  within  which  the  signal  is 
always  received  and  beyond  which  it  is  never  received. 
Instead,  the  situation  is  better  depicted  by  Fig.  3 
which  shows  theoretical  and  experimental  results  for 
the  correlation  between  the  signals  detected  by  two 
different  receivers  as  a  function  of  the  distance 
between  them.  When  the  unintended  receiver  is  moved  in 
a  direction  perpendicular  to  the  transmission  path,  the 


correlation  falls  off  more  rapidly  than  when  the 
unintended  receiver  is  moved  along  the  propagation 
path.  [9] 

While  the  theoretical  along-the-path  correlation 
falls  off  at  an  alarmingly  slow  rate,  the  experimental 
results  illustrate  two  key  points:  the  theoretical 
performance  is  difficult  to  achieve  in  practice,  and 
the  correlation  decreases  at  higher  operating  frequen¬ 
cies.  In  addition,  it  has  been  shown  [29]  that  fading 
reflections  are  responsible  for  a  significant  part  of 
the  correlation  at  large  separations  (trails  that  have 
been  distorted  by  winds  exhibit  broader  reradiation 
patterns  than  newly  formed  trails) .  By  communicating 
only  during  the  nonfading  portion  of  each  burst,  the 
size  of  a  footprint  can  be  kept  to  a  minimum. 

LPI  is  best  achieved  by  operating  in  a  burst  mode 
with  a  low  duty  cycle.  While  interception  by  a 
beyond-the -horizon  transmitter  is  then  governed  by  the 
probabilities  shown  in  Fig.  3,  a  meteor  burst  system  is 
vulnerable  to  line-of-sight  interception  unless  spread 
spectrum  techniques  are  employed. 

As  for  AJ  considerations,  line-of-sight  jamming  can 
be  easily  accomplished  in  view  of  the  losses  incurred 
by  the  desired  signal  in  the  course  of  being  reflected 


from  a  meteor  trail.  In  the  case  of  beyond-the-horizon 
j aiming,  the  j earner  can  distrupt  only  a  small  fraction 
of  the  traffic  if  he  relies  on  reflections  from  meteor 
trails  to  propagate  his  jamming  signals.  A  much  better 
stategy  is  to  use  extremely  high  power  (perhaps  in  a 
pulsed  model)  and  rely  on  ionoscatter  propagation. 
This  approach  can  be  defeated  by  operating  at  the 
higher  meteor  burst  frequencies,  since  ionoscatter 
falls  off  much  faster  with  increasing  frequency  than 
does  meteor  trails  reflections.  Finally,  a  jammmer  can 
attempt  to  spoof  the  system  by  transmitting  a  simulated 
probe  signal,  but  spoofing  can  be  defeated  by  judicious 
system  design. 


1.5.3  System  Perforfmance 

The  performance  parameters  of  greatest  interest  in 
commercial  applications  are  the  average  throughput  (in 
bit  per  second) ,  the  average  waiting  time  (to  transmit 
a  messge  of  specific  length) ,  and  the  probability  of 
error.  These  parameters  will  be  influenced  to  varying 
degrees  by  the  environmental  variables:  background 
noise  level,  external  interference,  geographic  loca¬ 
tion,  season,  and  time  of  day.  The  system  designer,  on 


the  other  hand,  can  exert  control  over  system 
performance  by  manipulating  his  design  variables  - 
transmitted  power,  RF  frequency,  instantaneous  data 
rate,  modulation  technique  and  antenna  gain. 

The  throughput  of  a  meteor  burst  system  depends  on 
factors  such  as  transmitted  power,  operating  frequency, 
communications  range,  received  noise  and  instantaneous 
data  rate.  The  increase  in  throughput  that  can 
ultimately  be  attained  is  limited  by  intersymbol 
interference  due  to  multipath  effects,  practical  bit 
synchronization  problems,  and,  of  course,  the  available 
transmitter  power  and  antenna  gains. 

In  addition  to  the  factors  under  the  control  of  the 
system  designer,  the  throughput  also  depends  on  such 
environmental  variables  as  time  of  day  and  season  of 
the  year.  These  variations  typically  involve  a  factor 
of  five  diumally  and  three  seasonally. 


1.6  Past  Systems 


A  landmark  meteor  trail  communication  system  was 
operational  in  the  1950s — the  Canadian  JANET  system 
(after  Janus,  the  Roman  god  who  looks  both  ways  at  once) . 


It  was  used  for  teletype  communications  between  Toronto 
and  Port  Arthur,  a  distance  of  approximately  1000  km 
[1] , [4] .  A  double  sideband  AM  duplex  channel  was  used  to 
transmit  and  receive  cosine-squared  shaped  pulses  in 
pulse  position  modulation  (ppm)  at  a  rate  of  650  bits  per 
second  (bps) .  Punched  paper  tape,  magnetic  tape,  and  a 
toriodal  magnetic  matrix  storage  core  were  used  as 
storage  media  to  compress  messages  in  time  for 
transmission  and  to  store  and  expand  the  messages  at  the 
receiving  station.  Despite  its  relative  antiquity,  the 
JANET  system  embodied  many  of  the  features  of  today's 
systems:  stored  digital  data  which  is  transmitted/re¬ 
ceived  in  bursts  after  carrier  detection,  duplex 
operation  at  VHF  frequencies  around  50  MHz,  and  low 
effective  duty  cycles  (about  0.1  for  the  JANET  system). 
Other  links,  such  as  a  one  way  Bozeman,  MONTANA-STanford, 
California  meteor  burst  link  were  also  operational  during 
the  1950s,  but  the  JANET  system  represents  the  first 
mature,  complete,  practical  hardware  implementation  of  a 
meteor  burst  system. 

Subsequent  experimental  tests  were  conducted  during 
the  1960s  and  1970s.  Notable  among  these  efforts  is  the 
experimental  work  summarized  by  Sugar  [4]  and  the 
implementation  of  the  COMET  (Communication  by  MEteor 
Trails)  system  which  has  operated  between  the  Netherlands 


and  Southern  France.  COMET  made  use  of  frequency 
diversity  and  Automatic  Repeat  for  reQuest  (ARQ)  which 
allowed  for  signal  repetition  in  the  event  of  no  radio 
path.  The  COMET  system  used  a  2000  baud  signaling  rate 
and  FSK  with  a  6  KHz  deviation.  The  storage  devices  were 
similar  to  those  of  the  JANET  system.  COMET  demonstrated 
the  practicality  of  meteor  burst  communications  under  a 
variety  of  conditions.  Typical  results  indicated  suc¬ 
cessful  transmission  of  a  50  baud  (average)  message,  150 
characters  long  with  a  delay  of  about  a  minute.  Worst 
case  delays  of  3  or  4  minutes  were  typical  [5].  Studies 
of  COMET-type  systems  have  continued  into  the  1980s 


2  MATHEMATICAL  MODEL  AND  CONCEPTUAL  FRAMEWORK 


In  this  chapter  the  basic  equations  that  model  the 
meteor  burst  channel  are  presented.  In  addition,  basic 
assumptions,  constraints  and  empirical  factors  are  stated. 


2.1  Power  Received  and  Related  Parameters 


The  power  received  from  the  trail  model  used  here  is 
from  the  classic  paper  by  Eshelman  and  Manning  [6]  and 
Hines  and  Forsyth  [16]. 


2.1.1  Power  Received  from  an  Underdense  Trail 


Meteor  trails  with  fewer  then  1014  electrons  per 
meter  of  length  are  referred  to  as  underdense. 
Eshelman  has  shown  that  a  simple  Fresnel  description  of 
the  process  is  normally  sufficient.  The  scattering 
then  is  specular  for  the  majority  of  recevied  signals 
that  is  the  incident  and  scattered  rays  make  equal 
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angles  6  with  the  axis  of  the  trail  and  the  received 
signal  coses,  in  effect,  from  the  principal  Fresnel 
zone  alone.  See  Fig.  4* 

The  principal  Fresnel  zone  is  centered  on  a  point  P 
in  the  trail  at  a  distance  Ry  and  Rr  from  the 
transmitter  T  and  receiver  R  respectively  such  that  R? 
+  Rr  is  a  minimum  for  the  trail  in  question.  The 
received  power  rises  rapidly  in  the  fraction  of  a 
second  taken  by  the  meteor  to  traverse  the  principal 
zone  and  reaches  a  peak  volume  at  time  t-=0 . 

Although  the  meteor  trail  is  formed  as  a  narrow 
column,  a  few  centimeters  in  diameter,  it  immediately 
begins  to  expand  by  diffusion.  As  the  diameter 
increases,  the  scattered  signal  suffers  increasingly 
from  a  destructive  interference  of  the  fields  scattered 
by  individual  electrons.  In  ideal  circumstances, 
neglecting  Fresnel  ripples  and  initial  diameter,  the 
received  power  decays  exponentially  in  time,  see  Fig. 
7.  The  above  is  summarized  in  the  next  equation. 
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L  1.1  C»«’A.l 


/,TCTCJ,r2sina(a)A.3exp^-^^j 

m  - - - -  2 


16 n2  RTRk{RT  +  R*){\  -  cos2 R sin 2 f ) 


~T T\«  exP 


32  n2P 
A.2sec2^ 


where, 


(2. 1.1-1) 


Pip  Transmitter  power . 

G<p,Gp  Transmitter  and  receiver  antenna  gains,  respec¬ 
tively. 

re  Classical  radius  of  an  electron  (2.82  10"15  m) . 
a  Angle  between  the  electron  field  vector  E  and 


q  Electron  line  density  of  the  trail. 
a.  Carrier  wavelength. 

ro  Initial  radius  of  the  trail  (.65  m.). 

*  Angle  of  incidence/reflection  of  the  transmitted 
plane  wave. 

Rp  Distance  from  transmitter  to  trail. 

Rr  Distance  from  receiver  to  trail. 

D  Diffusion  coefficient  of  the  atmosphere. 
fi  Angle  between  the  trail  and  the  propagation  plane 
formed  by  R«p  and  Rp. 


Figures  4,  5,  and  6  depict  the  geometrical  parameters. 


A  minimal  set  of  assumptions  is  applied  in  the 
literature  [7,11,21]  for  analysis: 
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1 .  )  The  trail  occurs  at  midpoint  between  receiver  and 
transmitter  i.e.f  • 

2. )  The  burst  occurs  at  an  altitude  of  100  km. 

3. )  The  trail  is  travelling  at  a  plane  perpendicular  to 
the  plane  formed  by  Rip  and  Rr  i.e.  p-x/2  .  in  addition, 
we  set  a- a/2. 


Under  these  assumptions  and  after  substituting  for 
the  physical  constants  we  rearrange  the  above  equation 
to  arrive  at  the  power  received  equation  for  the 
underdense  case. 


P(t)-Cuq7e-/B 


(2.1  .1  -2) 


where. 


Cv  -  2.5179581 •  10 
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D  -  -  ,  3  ,  (  33.359  \ 
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Rt 


H  -3.166-  10 


-3  1  sec  0 
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It  should  be  emphasized  that: 

1. )  P(0) ,  the  maximum  power  from  the  underdense  trail 
is  proportional  to  q2  ;  for  a  specific  trail  q  is 
fixed.  (The  random  character  of  q  is  discussed  later) . 

2. )  the  constant  of  proportionality,  incorporates 

the  effects  of  link  geometry,  transmitter  power, 


antenna  gains  and  carrier  wavelength.  For  a  given 
communication  system  Cy  is  constant  since  link 
distance,  transmitter  power,  antenna  gains  and  carrier 
frequency  are  set. 

3 . )  The  decay  time  constant  B  given  above  is  found  to 
be  randomly  varying  from  trail  to  trail.  The 
expression  in  the  classic  model  corresponding  to  a 
decay  time  constant  is  taken  as  the  average  decay  time 
constant  in  accordance  with  observations.  To  date, 
however,  researchers  have  assumed  all  underdense  trails 
to  share  the  same  decay  time  constant  to  f ascii itate 
their  analysis.  It  is  in  this  added  complexity  to  the 
model  that  our  treatment  deviates  from  published 
analysis  and  enhances  its  accuracy. 


4.)  The  decay  time  constant  B  and  the  electron  line 
density  q  are  functionally  independent  [1]  [7]  and 

hence  statistically  independent. 


5.)  For  a  given  communication  system  (i.e.  Cy  fixed) 
knowledge  of  the  electron  line  density  ,q,  and  decay 
time  constant  ,B,  completely  specify  the  underdense 


trail  behavior  in  time. 


Based  on  recent  empirical  data  [33]  the  pdf  for  q 


is: 


/,(<7)-Qq‘P  qmia<q<qu-  10Me/m  (2. 1.2-2) 

where, 

«? -(p-i 

p  -  1.6 

qmiB  -  Minimum  q  as  explained  below. 


*1  min  the  minimum  electron  line  density  q  a 
trail  must  posses  to  be  'seen'  by  the  communication 
system  i.e.  Cy  q2min  the  minimum  power  level  Pmin 
detected  by  the  communication  system.  In  the  context  of 
constant  transmission  rate  systems  (  as  disscussed 
later  )  PBin  is  the  power  level  that  corresponds  to  the 
maximum  allowed  bit  error  rate.  See  Fig.  8 
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2.1.3 


We  start  by  statin?  that  the  expected  value  of  the 


decay  tine  constant  B  is  veil  approximated  by: 


£(/?)»Z?-3.166I0~3-  ——^C—  —  sec  (2. 1.3-1) 


as  defined  in  section  2.1.1, 


The  random  decay  time  constant  B  is  assumed  to  be 


exponential  or  Rayleigh.  Both  assumptions  have  been 


shown  to  be  consistent  with  recent  experimental 


observation  [33].  In  general  the  Rayleigh  distribution 


will  be  used  since  for  very  small  B  (  extremely  fast 


decay  )  the  probability  is  very  near  zero.  This  is 


appealing  from  a  physical  viewpoint  and  consistent  with 


observation.  It  should  also  be  noted  that  even  if  such 


exceedingly  fast  decaying  trails  occur  they  are  of  no 


utility  for  communication. 


Exponential: 


(2.1  .3-2) 


where  9  is  defined  above. 
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Rayleigh: 


/.(b)--^e‘6i/2°2  (2. 1.3-3) 

aJ 

where, 


2.1.4  Joint  PDF  of  a  and  B  -  Pnderdense  Case 


Since  the  electron  line  density  and  the  decay  time 
constant  are  statistically  independent  the  joint  pdf  is 
given  by  their  product  i.e.  /,tf(g.b)-/,(g)-/,(b) 

where  the  individual  densities  are  given  in  the 
previous  sections. 

2.1.5  Power  received  from  Overdense  Trails 

Trails  with  more  10 14  electrons  per  meter  are 
often  termed  "overdense."  When  treating  them,  it  is  no 
longer  satisfactory  to  think  of  the  incident  wave  as 
penetrating  the  trail  without  serious  modification. 
Instead,  coupling  between  individual  electrons  plays  a 


prominent  part,  and  the  whole  of  the  scattering  process 
is  conceived  more  simply  in  terms  of  reflections  from  a 
cylindrical  surface. 

The  surface  in  question  surrounds  the  axis  of  the 
trail  at  a  radius  r,  where  the  electron  volume  density 
has  the  critical  value  normally  associated  with 
reflection  of  the  incident  wave  (cf.  reflection  from 
the  ionosphere) .  The  volume  density  within  the  trail 
exceeds  the  critical  value,  and  cannot  support 
propagation  in  the  usual  sense  at  the  pertinent 
frequency.  As  the  electrons  diffuse  outwards,  r 
increases  from  an  initially  small  value,  passes  through 
a  maximum,  and  falls  to  zero  as  the  volume  density  at 
the  axis  of  the  trail  decreases  below  the  critical 
value.  In  the  very  latest  stages  the  scattering 
reverts  to  the  underdense  type,  but  only  after  some 
time,  and  so  only  when  the  underdense  scattering  is 
extremely  small. 

The  scattering  process  is  again  specular,  the 
angles  of  incidence  and  reflection  being  equal.  The 
received  power  variation  is  given  by: 
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and  all  other  parameters  are  as  defined  for  the 
underdense  burst. 


Using  the  same  assumptions  as  for  underdense  case  ( 
a‘fi’ti/2  and  RT‘R„  )  and  some  algebra  we  rewrite  the 
above  as: 


+  I  t<ag-*(2.1.5-2) 


where. 
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As  in  the  underdense  case  we  have  rearranged  the 
classic  equation  for  ease  of  handling.  Note  that  C 0 


contains  the  transmitted  power,  antenna  gains, 
diffusion  coefficient,  carrier  frequqncy  and  the 
effects  of  link  distance.  Fig.  7 

In  addition,  from  eq.  (2. 1.1-2)  we  relate  CD  and  Cq 
by 


1  .2574816  10  29  — •  e^2***2*) 

A. 


(2.1  .5-3) 


2.1.6  Electron  Line  Density  Statistics 
-  Overdense  Case 


Based  on  recent  empirical  data  [33]  the  pdf  for  q 

is: 


/,(<7)-Qq'p  Q  muto  -  Q  -  Q  ho  “  10l6e/m  (2. 1.6-2) 
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Q  -(p-  l  ) <7 mmo  ‘  ^  jfl-  i 

p  -  1  .5 

g„Uo  -  Minimum  q  as  explained  below. 

qmino  is  the  minimum  electron  line  density  g  an 
overdense  trail  must  posses  to  be  'seen'  by  the 
communication  system.  Note  that  gnino  >  9Lo  the 
(physically)  lowest  possible  electron  line  density  for 
an  overdense  trail.  g^-lO14  e/m.  We  define  another 
constant  which  we  refer  to  later  as: 
Q'-l/tl-tqmino/qLoJP"1] • 


2.2  Link  Geometry 


In  this  section  we  present  the  relationship  between 
link  distance,  L,  and  station  to  trail  distance,  Rr  (»rr) 
and  angle  of  incidence/reflection  #.  See  Fig.  6 

The  station  to  trail  distance  is  given  by 


J{h  +  L2/ZRty  +  L7/4 


(2.2-1) 
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And  the  incidence/reflection  angle  #  is  found  using 


sec  (4) m  1 + 


where, 


(2.2-2) 


trail  altitude  (100  km) . 
great  circle  distance  between  stations, 
distance  from  transmitter  to  trail, 
distance  from  receiver  to  trail, 
radius  of  the  earth,  6400  km. 


The  derivation  of  the  above  is  done  in  section.  4.1 
Since  it  is  often  convenient  to  express  L  and  the 
resulting  Rr  in  thousands  of  kilometers  for  ease  of 
computation  we  then  have: 


RT-RRmyf{7\  +  (.01953)Z.j)2  +  Z.2/4 


sec  -  1  + 


(0 .2  ■*•  (  .039)Z.2): 


(2.2-3) 
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tor, 


s.t.  for  L-1000  km.  we  use  L*1  to  quickly  yield  RipmRR 
514  km.  and  sec2  #  -  18.5 
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In  this  section  we  are  concerned  with  quantifying  the 


trail  arrival  process  in  terms  of  communication  system 


parameters.  Since  the  arrival  is  Poisson  we  need  to 


estimate  the  average  time  between  bursts  or  equivalently 


the  number  of  bursts  per  second.  Oetting,  [9]  summarized 


the  results  of  the  COMET  system  showing  an  average  time 


between  bursts  of  4  to  20  seconds. 


It  is  often  neccessary  to  estimate  the  number  of 


arrivals  of  bursts  per  second  with  an  electron  line 


density  of  q  between  some  qi  and  q2.  We  let  A*  (q)  be 
the  number  of  trails/s  with  an  electron  line  density 


between  q  and  q+dq  e/m.  Using  Weitzen  [10]  we  write 


A'(q) 
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where  t  is  the  idealized  beamwidth  of  the  antenna  and  all 


other  parameters  have  been  previously  defined.  The 


number  of  trails/s  with  q  in  a  given  range  [qi,q2]»  A,  is 
found  by  integrating  the  above  density  over  q  yielding 


the  following 


V  V  V* 
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19.1667-  4>sin4>[RRS): 


L  92  J 


(2.3-2) 


note  that  for  overdense  trails  (q^  >  1014)  we  scale  the 


constant  19.1667  in  the  above  Eq.  to  be  equal  to  3. 


In  addition,  care  must  be  exercised  in  setting  qmin  since 


this  is  an  empirical  relationship.  For  example,  for  a 


1000  km  link  with  antenna  angle  of  46°  q^n  must  be  above 
1q13  e/m  for  the  formula  above  to  yield  reasonable  rates. 


Finally,  we  mention  again  the  dependence  of  the 


number  of  bursts/s,  A,  on  other  factors  such  as  global 


location,  time  of  day,  and  season.  These  factors  should 


appropriately  scale  any  value  of  A  used  for  analysis. 


For  signal  frequencies  appropriate  to  the  MBC  the 
dominant  noise  sources  are  galactic  and  man-made  in 
origin.  Consequently,  the  received  noise  power  spectral 
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density,  N0,  is  a  function  of  the  galactic  noise  picked 
up  by  the  receiver  antenna  and  the  receiver  thermal 
noise. 


No  is  modeled  as  follows: 


N a“  kT, 


- 1 

o 

\  y 3  1 

L  l*  \ 

.15  J  J 

(2.4-1) 


where. 


k  -  Boltzmann  constant,  1  .3805  •  10'23 J/K 


T0 “  290 


N0 -  4 • 10 


LKm  power  loss  between  the  antenna  and  receiver 
F-  receiver  noise  figure 
k  -  wavelength  in  meters 


We  choose  typical  values  for  Lp  and  F  as  1.3  and  2.5, 
respectively  to  yield  the  following  model  for  Nq 


-21 


80 


isj 


2.3 


+  2.5 


W/Hz 


(2.4-2) 
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The  bit  error  rate,  BER,  denoted  here  as  p^,  is  in 
general  related  to  the  bit  energy  E^  and  power  spectral 
density  of  the  received  noise  Nq  as 


*b(0-  g 


*6(0 


(2.5-1) 


where  g()  is  determined  by  the  modulation  technique. 


Since  Efc(t)«  P(t)/R(t)  where  R(t)  the  bit  transmission 


rate  may  be  time  varying  or  fixed  we  have  from  Eq. 


(2.5-1) 


*6(0-g 


*(0 


N0R(t) 


(2.5-2) 


or  altemtively,  we  can  solve  for  the  power 


pan  -N0-Ran-g"{pban) 


(2.5-3) 


where 


g‘‘(  )  is  the  inverse  of  g(-) 


t 


For  BPSK  we  have: 


/\(0 


Pj  0 
N0Rit)J 


Pit)  -  AV*(0[erfc*'(2 


For  BFSK  we  have: 


3  ANALYSIS  AND  RESULTS 


This  chapter  presents  and  explains  our  approach  and 
reasoning  to  the  use  and  analysis  of  the  MB  channel.  The 
significant  analytical  expression  to  be  used  as  design  and 
analysis  tools  are  given. 


3.1  Introduction  and  Approach 


As  stated  earlier,  our  aim  in  the  analysis  is  to 
arrive  at  a  set  of  analytical  relationships  for  important 
communication  parameters  such  as  burst  duration  and 
throughput  using  as  accurate  a  model  as  possible.  These 
relationships  will  serve  as  analysis  and  design  tools 
where  the  effects  of  system  parameters  and  constraints 
(power,  geometry,  BER)  on  throughput  can  be  assessed. 
Our  goal  is  to  estimate  the  optimal  average  throughput 
for  a  given  communication  system. 


Two  types  of  communication  systems  are  considered: 


1)  A  system  with  a  constant  transmission  rate  and  a 
time  varying  bit  error  rate  which  is  constrained  not  to 
exceed  some  maximum  allowable  value. 

2)  A  system  employing  a  variable  transmission  rate 
which  mimics  the  time  varying  power  in  order  to  maintain 
a  practically  constant  bit  energy  and  BER  for  all 
transmitted  bits. 

For  the  constant  bit  rate  system,  since  the  power  is 
time  varying,  the  probability  of  bit  error  will  be  time 
varying  in  a  fashion  dictated  by  the  relevant  modulation 
function.  Specifically,  for  the  underdense  burst  we  have 
an  exponentially  decaying  power  in  time  and  consequently 
the  BER  is  monotonically  increasing.  Whether  underdense 
or  overdense  we  see  the  need  then  to  impose  a  ceiling  on 
the  BER.  For  such  systems,  where  the  probability  of  bit 
error  is  time  varying,  we  specify  a  maximum  allowable 
BER,  Pbrnax*  such  that  all  tranmitted  bits  shall  have  a 
BER  less  than  or  equal  to  Ppmax*  This  constraint  on  the 
probability  of  bit  error  for  constant  transmission  rate 
systems  implies  the  existence  of  a  minimum  power  level 
Pmin  such  that  all  transmitted  bits  have  a  received  power 
greater  than  or  equal  to  PBin.  Such  is  obvious  when 
remembering  that  the  modulation  function  relating  power 
to  BER  is  an  increasing  function  for  decreasing  argument. 


wwnrcvwuiTj : 


mm  w.  W77W 
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In  short,  Pfc(t)  ■  P^max  w^en  P(t)  «  pnin*  The  time  for 
which  the  power  received  from  the  meteor  burst  exceeds  a 
prescribed  threshold  (Pain)  i®  the  burst  duration  tg 

(also  duty  cycle  [1]). 

Clearly  the  duration  of  time  for  which  the  received 

power  P  is  greater  than  PBin  varies  from  one  meteor  burst 

to  the  next.  For  a  constant  bit  rate  system  the  burst 
duration  is  crucial  in  determining  the  number  of  bits 
transmitted  during  a  given  burst  and  ultimately  in 
estimating  throughput.  We  shall  therefore  analyse  its 

relationship  to  system  parameters,  derive  its  statistics 
and  compare  it  with  current  dogma. 

As  we  increase  our  set  constant  bit  rate  the 
corresponding  P^,,  has  to  be  increased  resulting  in 
shortening  the  effective  duration  of  transmission.  This 
clear  tradeoff  between  transmission  rate  and  duration  of 
transmission,  tg,  suggests  an  optimal  choice  for 
transmission  rate  (or  equivalently  an  optimal  tg)  so  as 
to  maximize  the  throughput.  The  link  between  the  number 
of  bits  transmitted  per  burst  and  throughput  is  done 
through  the  average  arrival  rate  statistics  and  trail 
shape  statistics  (variation  in  q  and/or  B)  .  The 
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judicious  choice  of  transmission  rate  for  the  given 
system  provides  us,  therefore,  with  the  optimal 
throughput . 

The  second  type  of  communication  system  considered 
here  is  a  variable  transmission  rate  system.  Under  such 
a  system  the  bit  rate  mimics  the  time  behavior  of  the 
received  power  (or  equivalently  stated  the  bit  duration 
varies  in  a  reciprocal  manner  to  the  received  power 
variation)  so  as  to  maintain  a  constant  bit  energy.  The 
direct  outcome  of  such  approach  is  a  constant  probability 
of  bit  error.  More  importantly,  we  note  that  we  are  not 
restricted  by  the  need  to  operate  above  a  prescribed 
threshold  of  received  power;  for  as  the  power  diminishes 
we  simply  compensate  by  increasing  the  bit  duration. 
This  adaptive  approach  implies,  theoretically  at  least, 
for  the  underdense  burst  an  infinitely  long  time  of 
transmission.  This  of  course,  is  offset  by  the  practical 
limitation  on  bit  duration  and  the  need  to  maintain  a 
fairly  constant  BER.  The  assumption  here  is  that  the 
hardware  is  fast  enough  to  vary  the  transmission  rate 
(bit  duration)  from  one  bit  to  the  next.  Since  the  power 
in  most  cases  diminishes  during  the  duration  of  the  bit, 
for  long  bit  duration  we  deviate  from  the  requirement  to 
maintain  a  constant  bit  energy  and  constant  BER.  It  is 
therefore  necessary  to  take  this  fact  into  account  in 


determining  the  duration  of  transmission.  This  criterion 
of  allowing  only  a  small  drop  in  the  assumed  constant  bit 
energy  imposes  a  restriction  on  the  maximum  bit  duration 
which  in  turn  limits  the  duration  of  transmission.  As  in 
the  previous  case  we  estimate  the  throughput  by  averaging 
over  the  ensamble  of  trails  and  incorporating  the  effects 
of  trail  arrival  statistics.  Finally,  we  express  the 
correct  packet  probability  in  terms  of  the  modulation 
function,  packet  length  and  communication  system 
parameters . 


3.1.1  Protocol 

Both  stations  have  receiving  and  transmitting 
capabilities.  The  station  which  acts  as  transmitter  of 
data  listens  continuously  for  a  continuous  tone  which 
is  being  sent  from  the  receiver.  Upon  detection  of 
this  tone  both  the  presence  and  the  strength  of  the 
channel  is  known  to  the  transmitter  which  immediately 
commences  transmission.  The  situation  is  completely 
symmetrical  with  respect  to  the  receiver.  Clearly,  the 
beginning  of  the  burst,  i.e.  availability  of  the 
channel,  is  known  to  the  transmitter  within  a  neglibile 


propagation  delay  (less  than  3.5  nsec  for  a  1000  km 
link) .  The  continuous  tone  used  for  monitoring  the 
channel  is  at  a  different  frequency  than  the  data  for 
ease  of  gating  and  synchronizing.  For  constant  bit 
rate  systems  the  essential  feature  of  this  scheme  is 
the  ability  of  the  transmitter  to  discern  the  beginning 
of  the  burst.  For  variable  bit  rate  systems  the 
strength  of  the  channel  as  gauged  by  the  transmitter  in 
addition  to  the  burst 1 s  start  are  important  for 
adaptive  bit  rate  determination. 

The  above  description  is  the  rudimentary  layer  upon 
which  additions  and  refinements  could  be  made. 


3.2  Improvements 

The  advantages  of  the  analysis  herein  stem  basically 
from  considering  a  more  complex  MBC  model  than  found  in 
the  literature.  Most  researchers  [9,  21]  while  acknowl¬ 
edging  the  significant  contribution  to  communication  of 
the  rare  yet  robust  overdense  burst  have  tended  to 
concentrate  on  the  underdense  burst  phenomenon  because 
its  simpler  exponential  form  is  more  mathematically 
tractable.  Able  [11] ,  however,  does  analyze  the 
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overdense  trail  by  a  piecewise  linear  (first  order) 
approximation  resulting  in  conservative  channel  duration 
estimation.  We  choose  a  second  order  approximation  to 
improve  the  fit  of  the  overdense  transcendental  power 
time  function.  This  yields  better  results  for  channel 
duration  for  the  overdense  case.  For  the  exponential 
case,  researchers  to  date  have  assumed  a  fixed  decay  time 
constant.  Hampton  [21]  in  his  analysis  of  an  MB 
broadcast  system  used  a  constant  decay  rate  which  was  the 
conditional  expected  time  constant  given  a  fixed  burst 
duration.  This  approach,  however,  implies  that  the  decay 
time  constant  is  dependent  upon  an  average  of  a  function 
of  the  electron  line  density  whereas  it  is  known  that 
the  time  constant  is  independent  of  the  electron  line 
density.  To  enhance  the  accuracy  of  our  model  we  chose 
to  adhere  to  empirical  results  and  assume  the  decay  time 
constant  to  be  a  random  quantity  independent  of  the 
electron  line  density. 

Some  researchers  [34]  have  further  simplified  the  MBC 
model  by  assuming  a  fixed  electron  line  density  for  all 
exponential  bursts.  This  is  contrary  to  experimental 
observation.  Here,  as  in  the  case  of  the  exponential 
decay  time  constant,  we  incorporate  this  random  parameter 
into  our  model. 


In  analyzing  throughput  we  compare  two  systems: 
constant  bit  rate  system  and  variable  bit  rate  system  and 
derive  the  average  throughput  for  both.  For  constant  bit 
rate  we  derive  the  expression  for  the  optimal  average 
throughput.  This  optimal  average  throughput  is  achieved 
by  finding  the  "best"  bit  rate  to  be  used  for  the 
communication  system  so  as  to  maximize  throughput  for  the 
underdense  burst.  This  is  not  to  be  confused  with  the 
optimal  bit  rate  for  a  given  burst  that  maximizes  the 
number  of  transmitted  bits  for  the  burst  under 
consideration . 

In  summary  we  note  that  the  bedrock  of  a 
communication  system  using  the  MBC  is  its  model.  To 
enhance  the  analysis  we  used  a  model  with  as  few 
simplifications  as  possible  while  maintaining  mathemati¬ 
cal  tractability. 

3.3  Sample  System 

For  quantitative  appreciation  of  the  results  we  need 
to  assume  some  values  for  a  communication  system.  The 
set  of  these  assumed  parameters  values  will  comprise  what 
will  heretofore  be  refered  to  as  our  sample  system.  This 
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sample  communication  system  represents  practical  or 
average  values.  Chapter  2  contains  the  definition  of 
the  parameters. 


1000  km  link  distance.  From  section  2.2 
ve  have  for  the  next  two  parameters: 


sec (*) 

rr 


=  18.5  and 


514  km. 


1000  v 


10  dB. 


Pbmax 


10  dB.  And  for  antenna  beamwidth  ve  use 


=  45.84°  (.8  rad) 


=  6  m.  (50  MHz) 

«  8  m2/s 

*  10”5  Max.  BER  for  constant  bit 
rate  system. 


Using  BPSK  modulation  ve  get: 


g~1(Pbmax)  *  9 


N0  -  4  •  10*21  80(^)  +  2.5  W/Hz 


See  discussion  above  Eq.  (2.4-2). 
4.89  10“20  at  50  MHZ. 


For  comparison  the  BER  for  variable  rate  system,  Pb» 


is  set  equal  to  Pbmax 

RminV  ■  1000  b/s  (minimum  variable  rate  allowed) 


3.4  Constant  Bit  Rate  and  Constrained  BER  System 


3.4.1  Channel  Duration 

Channel  duration  or  burst  duration,  tB  is  the  time 
for  which  the  power  received  from  the  burst  exceeds  a 
given  minimum  value  Pjnixi*  is  also  the  duration  of 

data  transmissions  for  a  given  burst  where  all  bits 
possess  a  probability  of  error  less  than  some  maximum 
allowable  BER,  Pbmax*  pmin  *s  determined  thru  the 
modulation  function  by  the  following  parameters: 
Transmission  rate  (non-time-varying) ,  power  spectral 
density  of  the  received  noise,  and  the  maximum 
allowable  bit  error  rate.  The  burst  duration  tB  is 
determined  by  Pmin  and  cu  which  incorporates  system 
parameters  such  as  transmitter  power,  link  geometry 
effects,  carrier  wavelength  etc. 
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3. 4. 1.1  Burst  Duration  Statistics  for 

Constant  Bit  Rate  --  Underdense 

The  channel  duration ,  -  burst  duration  -  is 

the  time  of  transmission  for  a  given  burst  i.e.  the 
usable  portion  of  the  burst,  assuming  a  constant 
transmission  rate,  R,  some  given  maximum  allowed  BER, 
and  power  spectral  density,  Nq  ,  see  Fig.  11. 

For  such  conditions  for  the  underdense  burst  we 
have: 


t „  -  B  •  In 


(3.4.1  .1  -  1) 


where 


9  min 


P 

1  mm 

C  u 


R  min  “  N0’  R’  g''{Pbnax) 


and 


Ctf  -  2.5179581  10 


-32 


.  Rt 


or  equivalently 


tB  -  B  •  In 


( 


C  u 


{N0-R-g-'{Pbmax) 


(3.4.1  .1-2) 


€ 1 


We  now  average  over  the  ensemble  of  underdense  bursts 
(averaging  over  q  and  b)  to  yield  the  average  burst 
duration: 


•}  T5  l 

fi--— ^•^•(l-x(l-lnx)]  (3. 4. 1.1-3) 


where 


# 

\  Qu  )  _  Cvql 


5-3.166  10~3  •  ^  S^C  *  qv-  10M  p-1.6 


ft  max  Q U 


C  uQu 


No9  '{P  bmax  } 


For  a  given  communication  system  Cy,  R,  Nq,  Pbmax, 
5  and  g( •)  are  specified.  R^  is  the  normalized  (  with 


respect  to  the  highest  possible  rate  for  an 


underdense  burst)  bit  rate. 


I . **#.**#  ’  •«#  Ul  l.»l».yi-l^Ll'll.tli-tlt.**»-<,*.a<4jilt«  t  «»> 


,1a>  Id 
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Fig.  14  plots  the  above  the  above  for  our  sample 
system  (Sec.  3.3)  as  a  function  of  %.  (RN«l — R-86.5 
kb/s) 


Ikg.  variance  of  tg  is  given  bv: 


o;  - 


1 


(P-1)2  (1  -x): 


[c,a2  +  c2Ji2]  (3. 4. 1.1-4) 


where 

ci  -  2  -  x (3  +  y2)  +  x2(l+y2) 
C2  *  1  -  x(2+ln2(x))  +  x2 
y  *  1  -  ln(x) 

x  from  previous  equation 


and 


o,_  - 


'*  (p-l)(l-x) 


- — -^kc,  +  c  j  •  B 


(3.4.1  .1  -5) 


For  B  Rayleigh 


A:  -  -  -  1  ,  For  B  exponential  k  -  1 


The  nth  moment  of  ts  is  given  bv: 


p  - 

•'(p-ir’i-x 


l-xi(-i)1 


i-0 


»n‘(x)j 
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where 

x  and  p  are  from  above  eq.  (3. 4. 1.1-3) 

P'  “(P  “l)/2 

*>min  »  tx  /  [  (l/p«)  ln(l/x)  J 
Q'  -  l/(l-x) 

For  the  case  of  x«l  i.e.  Pmin  «  Cg  , 

essentially  ignoring  Pbaax  (considering  all  physical 
trails)  we  have  for  the  solution  of  the  above 
integrals  (  bj,in  -  0  ) : 

B  assumed  exponential: 

/«,(*, )-2y-*0(z)  (3. 4. 1.1-8) 


... ... 


'«  $  '4  I* 


'jr^mwrvrv1.  w  vjw 


and 


*■..(«,)- l-zX.(z) 


wherez  - 
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t ,  >  0  (3. 4. 1.1-9) 


where,  K  q  and  are  the  modified  Hankel 

functions  of  order  0  and  1  respectively. 


B  assumed  Rayleigh: 


F,(tl)-l-2  bexp(-b3-t,/b)db  (3.4.1.1-10) 

Jo 

-1-2/, (I,) 

where  fj.(ti)  is  a  tabulated  function  [35].  See 
Fig.  9. 

The  above  results  are  presented  in  greater  detail  in 
ref.  [12] 


3. 4. 1.2  Burst  Duration  Statistics  for 

Constant  Bit  Rate  —  Overdense 

An  exact  solution  to  the  roots  of  the  power 


received  equation  -  overdense  -  is  not  known.  The 
difference  between  these  time  instances  yields  the 
duration  of  the  channel.  We  can  approximate  the 
power  function  in  a  piecewise  linear  fashion 
(comprised  of  two  straight  lines)  or  alternatively 
use  a  second  order  approximation.  The  second  order 
approximation  used  here  improves  the  accuracy  of  the 
results . 

The  results  are  stated  without  the  lengthy 
derivation.  For  derivation  see  section  4.3 


tB=.95aqyJ  1  qn,no$q<qU0-  1016e/m  (3.4.1  .2-  1) 

where 

a  and  Cq  are  from  Eq  (2. 1.5-2) 

<3T  *  (e/a)  (Pttin/co  )2  =  <5U0  RN02 

rN0  *  R/Rmaxo  »  normalized  bit  rate. 

Rmaxo  “  C0[a  quo/e]  ’  VNo9-1  (Pbmax) 

maximum  possible  bit  rate  for  any 
overdense  burst. 

qmino  “  ®ax  (<JT/<lLO>  qLO-1014  quo  “  1016 

e/m. 

e  is  the  natural  logarithm  base. 
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Note  that  qjx)  is  the  minimum  (physical)  electron  line 
density  that  can  support  an  overdense  burst.  qT  is 
specified  by  our  choice  of  system  parameters  and  bit 
rate  and  may  not  correspond  to  a  real  level  of 
overdense  electron  line  density.  See  discussion 
following  Eq.  (3. 4. 1.2-3).  See  Fig.  12.  but  note 
that  in  the  figure  qi>«lLo* 

We  can  also  relate  qT  to  q^^  recalling  the 
definition  of  Pmj.n  -  Cy  q2min  ; 


Qt  - 


(3.4.1 .2-2) 


where  the  parameters  Cv,  C0,  a  and  qmin  are  defined 
in  Eqs.  (2. 1.1-2)  and  (2. 1.5-2)  and  (3. 4. 1-1) 
respectively. 

Average  Burst  Duration  -  Overdense: 

The  average  burst  duration  for  the  overdense  case  is 
found  by  averaging  over  the  ensemble  of  bursts, 
yielding 


For  the  case  where  RjfO  <  .1  (q^  >  q^)  we  have 


v  v  v  v"*  v  v  v  v  h."  v  v  v  v"  v  k_  '-’  "*"  v  w^  t” 


TFTT*X 
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Fig.  15  plots  the  above  for  our  sample  system  (Sec. 

3.3)  as  a  function  of  RN0.  (Rno“°*  1— R-256  kb/s) 

For  many  practical  systems  ve  generally  have 

qT<<qLo<<quo  i«®.  <3mino  “  ^L0»  Q"(9L0)l^2/2 
resulting  in: 

? t  “  •  95 a  •  •J Q uoQ lo  (3. 4. 1.2  —  4) 

and  for  q<j  <  .01  quo  v®  have  for  the  burst  duration 
variance : 

*?,- -95y/g^g^a- .9Syl031  (3. 4. 1.2-5)  ! 

For  our  sample  system  definsd  in  Sec.  3.3  the 
standard  deviation  of  tB  «  11  sec. 

I 

Burst  Duration  pdf  and  cdf  -  Overdense: 


f  ..(«,)■ 


tB: 

« 

1  - 

^  2 

I 

\ 


p-  1 


9r[l-Vl*(«,/r)a]J 


^  I  min  —  ^  I  —  ^  I  max 


(3.4.1  .2-6) 


: 


swww»rww  pi  iii  i.«j — mimij* 


where 


p«1.5  and 


‘95aqmUoJ  1  ~  ~~  qmino  m  max {g7 , q L0) 
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(3.4.1  .2-7) 


where 


and  all  other  parameters  are  defined  by  previous 
equations . 


.4.2  Bits  per  Given  Burst 


Bits  per  given  burst  under  constant  bit  rate  and 


constrained  BER  is  the  number  of  bits  transmitted  for 
the  burst  whose  q  and/or  B  are  known.  Since  the  bit 
rate,  R,  is  assumed  constant  the  number  of  bits,  NB,  is 
simply  the  product  of  R  times  the  burst  duration  tB  as 
defined  in  previous  sections.  Clearly,  the  level  of 
power  Pmin  which  reciprocally  affects  tB  is  dictated  by 
the  modulation  function  and  in  general  increases  with 
increasing  bit  rate  R.  As  a  result  we  see  the 
trade-off  between  the  bit  rate,  R,  and  transmission 
duration,  tB,  such  that  increasing  one  decreases  the 
other.  This  suggests  the  existence  of  an  optimal  bit 
rate,  R*B,  (or  equivalently,  optimal  transmission 
duration)  that  maximizes  the  number  of  bits  transmitted 
during  the  given  burst.  Clearly,  using  any  bit  rate 
other  than  the  optimal  R*B  will  result  in  under-uti¬ 
lization  of  the  burst.  Alternatively  stated,  given  a 
bit  rate  one  should  transmit  only  for  as  long  as  the 
optimal  burst  duration  dictates.  The  optimal  bit  rate 
for  the  given  burst  is  then  used  to  calculate  the 
maximum  number  of  bits,  NB*  ,  that  can  be  transmitted 
during  this  burst.  A  ratio  is  then  obtained  for  the 
same  burst  between  the  number  of  bits  transmitted  using 
variable  bit  rate  (BER  constant  -  Pbmax)  and  the  NB* 
using  constant  bit  rate.  This  ratio  has  been  shown  by 
Abel  [11]  to  be  around  2.5.  Using  the  "best"  (highest 
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q)  burst  for  overdense  and  underdense  trails  bounds  can 
be  found  for  the  maximum  possible  Nb*  (constant  rate 
system)  and  maximum  transmitted  bits  using  variable 
rate. 

Two  points  however,  must  be  kept  in  focus: 

1)  Since  R*  varies  from  trail  to  trail  and  a  priori 
knowledge  of  the  trail  behavior  is  impossible  to 
obtain,  the  optimal  (constant)  bit  rate  for  the  system, 
i.e.  for  all  possible  bursts,  is  still  unknown. 
Par anthet ica  1  ly ,  we  note  that  attempting  to  find  the 
burst's  shape  from  knowledge  of  its  beginning  is  both 
wasteful  of  valuable  transmission  time  and  unreliable 
since  many  bursts  deviate  from  ideal  behavior.  As  will 
be  shown,  the  important  parameter  for  throughput 
determination  is  not  Nb  —  the  number  of  bits  for  a 
given  burst  —  but  rather  the  number  of  bits  per  burst 
averaged  over  the  ensemble  of  all  bursts  since  the 
latter  is  directly  proportional  to  the  average 
throughput . 

2)  While  it  is  true  that  more  bits  can  be  transmitted 
using  variable  bit  rate  during  the  burst  rather  than 
employing  constant  bit  rate,  the  situation  is  somewhat 
analogous  to  an  apples/oranges  comparison.  For  the 
underdense  burst  and  constant  bit  rate  all  bits  except 
the  last  one  will  possess  a  BER  smaller  than  Pbmax  — 


the  worst  BER,  whereas  for  variable  bit  rate  all  bits 
will  possess  the  same  BER  -  P^max*  Consequently,  the 
probabiliiity  of  having  an  error-free  packet  of  data  is 
higher  for  constant  bit  rate  than  for  variable  bit 
rate. 

We  present  the  analysis  for  the  given  underdense 
burst  by  first  considering  constant  bit  rate  with 
constrained  BER  and  then  considering  variable  bit  rate 
system  with  constant  BER. 

Constant  Bit  rate  system 

Given  an  underdense  burst  -  q  and  B  are  assumed  known. 
Pbmax  is  specified  for  the  system  and  R  is  the  constant 
bit  rate.  Let  N3  be  the  number  of  bits  transmitted  on 
this  burst.  All  bits  have  BER<P]Mnax.  All  bits  have 
power  P> Pmin*  From  the  discussion  on  burst  duration 
and  the  above  definition  of  N3  we  have: 

N„-R-t ,  (3. 4. 2-1) 

Using  eq.  (3. 4. 1.1-2)  for  tB  we  rewrite  the  above  as 

N R-  B-\n[R0/R)  (3. 4. 2-2) 

where  Rq  *  Cu  q2  /  Nq  g"1 (Pbmax) 


Note  that  R0  is  a  constant  set  by  the  choice  of  system 
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parameters  and  physical  characteristics  (1)  of  the 
given  burst.  From  the  above  it  is  clear  that  for  a 
given  burst  a  different  choice  of  bit  rate  R  yields 
different  number  of  transmitted  bits,  Ng.  To  find  the 
bit  rate  Rg  that  maximizes  NB  we  simply  set  its 
derivative  with  respect  to  R  to  zero  and  solve 

-  0  -  \n(R0/ R.)  -  1  (3. 4. 2-3) 

dR 


or  R*g  =  Ro/e  where  e  is  the  natural  logarithm  base. 
Using  the  last  result  in  eq  (3. 4. 2 -2)  yield  Ng*  the 
maximum  number  of  bits  for  this  burst. 


(3. 4. 2-4) 


where  Rq  **  Cu  q2  /  N0  g-1(Pbmax) 


From  the  above  we  see  that  P*min  the  minimum  power 
corresponding  to  R*g  is  equal  to  P(0)/e.  This  implies 
that  the  burst  duration  (transmission  time)  i.e.  the 
time  it  takes  to  drop  to  P*min  from  P(0)  is  equal  to 
one  decay  time  constant  B.  Equivalently,  then,  the 
above  analysis  states  that  for  a  given  bit  rate  R  the 
optimal  time  of  transmission  T*g  must  be  set  to  B  to 
yield  the  above  maximum  transmitted  bit  N*g.  This 
result  has  been  demonstrated  by  Abies'  work  [11]  as 


ell. 


Variable  £it  B.gt9  gystem 


For  the  given  underdense  burst  we  have  Pb(t)“Pbmax 
for  all  bits.  Under  such  conditions  the  variable  bit 
rate  R(t)  is  given  by 


R(t)  *  Roe-t/B  (3. 4. 2-5) 
where  R0  -  C^/g"1  (Pbmax) 

Note  that  the  power  to  transmission  rate  ratio  which 
specifies  the  bit  energy  is  constant  and  hence  the  BER 
is  maintained  constant. 

Since  the  burst  is  assumed  to  be  known  a-priori  we 
can  compute  the  bit  duration  (inverse  of  bit  rate)  for 
each  bit  exactly  such  that  the  integral  of  power 
received  over  the  individual  bit  duration  i.e.  the  bit 
energy  is  the  same  for  all  transmitted  bits.  See  Fig. 
13. 


The  number  of  bits  for  variable  rate  is  given  by 


•jCTJC. 


(3. 4. 2-6) 


Ntv-  R(t)dt-  R0- B 
Jo 

Comparison  ol  TK2  Systems 

We  now  compare  constant  bit  rate  system  using  the 
optimal  R  (with  BER  <  Pbmax)  for  the  particular  burst 
versus  variable  bit  rate  which  mimics  the  bursts  power 
behavior  (such  that  all  bits  possess  the  same  BER, 
Pbmax)  •  For  each  of  the  above  we  found  the  number  of 
bits  transmitted  N*B  and  NBV  respectively.  Taking  the 
ratio: 


NBV/N*B  -  e 

We  see  an  improvement  of  2.7  for  the  utilization  of 
variable  bit  rate  over  constant  bit  rate.  We  stress 
however  that  the  necessary  priori  knowledge  (or 
estimation)  of  burst  characteristics  is  unreliable.  In 
addition,  although  variable  bit  rate  results  in  a 
higher  number  of  transmitted  bits  for  a  given  burst 
versus  constant  bit  rate  there  is  a  price;  all  bits 
have  a  BER  -  Pbmax  f°r  variable  bit  rate  system  whereas 
for  the  constant  bit  rate  system  Pbmax  is  a  worst  case 


(last  bit  BER) .  This  fact  portends  worse  probability 
of  correct  packet  transmission  for  variable  rate  system 
than  constant  rate  system  tinder  the  above  conditions. 


3.4.3  Throughout  for  Constant  Bit  Rate 

We  are  given  a  communication  system  with  a 
specified  maximum  allowed  BER  Pfemax'  constant  transmis¬ 
sion  rate  R,  received  noise  power  spectral  density  Nq 
and  general  system  parameters  (link  distance,  frequen¬ 
cy,  etc.)  as  defined  by  Cy  or  Cy. 


To  express  the  throughput  T  of  the  system  we 
introduce  Nyi  as  the  number  of  bits  transmitted  on  the 
i'tb  burst  and  m,  as  the  number  of  bursts  occurring  in 
the  period  of  r  seconds.  Since  the  bursts  are 
independent  of  each  other,  and  assuming  no  overlap  of 
bursts,  we  have: 


T  (3. 4. .3-1) 

‘  i- 1 

-  A-N, 


-**1 


where  A  is  the  average  number  of  bursts/s  and  N,  is  the 
expected  number  of  bits  per  burst  as  averaged  over  the 
ensemble  of  all  bursts'  profiles.  (The  subscripts  U  or 
0  are  added  to  distinguish  between  underdense  and 
overdense  cases . ) 

For  underdense: 

Nbu~J  j  N  q(q)f  B(b)dqdb  (3. 4. 3-2) 

where  fx  ()  is  the  pdf  of  x. 


For  overdense: 

/ N,(g)/,(<7)d<7 


(3. 4. 3-3) 


Note  that  the  above  three  equations  are  dependent  on 
the  choice  of  bit  rate  R. 


3. 4. 3.1  Optimal  Average  Throughout  for  Constant: 

Bit  Rate  System  —  Underdense. 
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In  this  section  the  optimal  average  system 
throughput  for  the  underdense  burst  is  given.  This 
is  the  maximum  average  number  of  bits  per  unit  time 


for  the  *best'  transmission  rate  R  (and  other  system 
constants)  assuming  only  underdense  bursts  are 
present. 

This  optimal  throughput  is  found  by  first  deriving 
the  dependence  of  the  average  throughput  (Eq. 
3. 4. 3-1)  on  the  bit  rate  R,  and  solving  for  the 
optimal  bit  rate  R*  that  would  yield  the  desired 
optimal  throughput  T*.  The  solution  is  done  with 
respect  to  the  normalized  bit  rate,  RN,  which  yields 
a  result  independent  of  a  particular  choice  of  system 
values . 


For  the  given  system  Pbmax*  N0,  R  and  cu  are 
specified.  We  have  defined  RN  the  normalized  bit 
rate  as: 


R 

R  max 


Q  m 
<7 


n 


2 
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where 


R 


max 


£_  uQu 

No  Q  '  (  P  b  max  ) 


(3.4.3. 1  -  1) 


the  maximum  bit  rate  achievable  in  a  given  system 
using  underdense  bursts. 


Omitting  the  cumbersome  derivation  we  present 
the  following  results:  (for  the  relevant  derivation 
see  Sec.  4.2)  * 

The  average  number  of  bits  per  underdense  burst 
Ntu-R-T,  can  be  expressed  in  terms  of  the  normalized 
bit  rate  using  Eg.  (3.4. 1.1-3)  as: 


A7  „„  - 


P-1  ““1-K& 


z[RN-RpN"  +  RpN"-ln{RpN)]  (3.4.3. 1-2) 


where 


p'«(p  -  l)/2  ,  p*l. 6  and  Rmax  *s  defined  above. 


Note  that  in  the  above  expression  the  bursts 
considered  were  those  with  an  electron  line  density  q 
between  q^in  (as  defined  by  Pmin)  and  qu  (the 
physical  limit) .  The  arrival  rate  for  such  condition 
is  given  by  Eq.  (2.3-2)  where  qx  is  set  to  qmin  and 
<12  is  <Ju*  Together  with  Eq.  (3. 4. 3. 1-1)  we  have  for 
the  arrival  rate  of  underdense  bursts: 
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where, 


19.1667-*sin*(/?,£)2 

/I  - - 

Qu 


Qu 


10 


14 


(3.4.3. 1  -3) 


The  expression  for  throughput  using  underdense  bursts 
is  given  by  the  product  of  the  last  two  Egs.  as 
specified  by  Eq.  (3. 4. 3-1).  The  result  is  a  function 
of  the  normalized  bit  rate  which  is  just  a  scaled  bit 
rate.  In  the  product  we  substitute  for  Rmax  from  Eq. 
(3. 4. 3. 1-1)  and  for  Cy  and  B  from  Eq.  (2. 1.1-2).  For 
a  given  set  of  system  parameters  we  thus  have  the 
throughput  (using  underdense  bursts)  for  any  choice 
of  (normalized)  bit  rate  as: 


Ty-  3.056-  10 


PtGtGrS2*  sin*sec2*expf-^-Vl; 

-  19  _  _ V  1‘wi'lj 

(p-  1  )N 0-  g'x[pbmax)-  D  ■  Rr 
•]  R  N  ~  R  N  I 


1  -RV 


'[l-RpN+p'RpN-\n[RN)] 


(3.4.3. 1  -4) 
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where 


P 

Rp  and  * 


PT 

Gp,GR 

A. 

D 

N0 


9”1(pbmax) 


is  an  empirical  constant  ,  p'«(p-l)/2 
(taken  here  as  1.6,  Section  2.1.2) 
idealized  antenna  beamwidth. 
geometric  parameters  detrained 
by  link  distance  (from  section  2.2) 
transmitter  power. 

transmitter  and  receiver  antenna  gains, 
carrier  wavelength. 

diffusion  coefficient  of  the  atmosphere, 
power  spectral  density  of  the  noise. 
(Sec.  2.4) 

ratio  of  minimum  bit  energy  to  N0. 


Note  that  the  equation  is  written  in  two  parts;  the 
first  incorporates  all  the  parameters  of  the 
particular  system  and  the  second  part  reflects  the 
variation  in  throughput  with  choice  of  bit  rate.  The 
above  throughput  (divided  by  about  3200  bits/s)  for 
our  sample  system  (defined  in  Sec.  3.3)  is  plotted  in 
Fig.  16.  Note  the  maximum  is  214  b/s. 

We  would  like  to  find  from  the  above  equation 
the  optimal  normalized  bit  rate  R^*  which  yields  the 


maximum  throughput  T^*  .  This  optimal  normalized  bit 
rate  for  any  given  system  is  found  numerically  to  be 
(see  Fig.  16) : 


Rn*  -  .07 

or  equivalently  the  optimal  bit  rate, 

R*  =  .07  Raax  (3. 4. 3. 1-5) 

resulting  in 


Ty -  3.433 • 10 


-20 


^rGTG*C20sin(0)sec2(0)expf-^^- V 

RTN0g-'[Pbmar)D 


(3.4.3. 1  -6) 

where  all  the  parameters  are  defined  above. 


The  above  expression  for  the  maximum  throughput  Ty* 
allows  us  to  quantify  the  affects  of  various  system 
pareuneters  and  constraints  on  the  best  throughput 
using  a  constant  bit  rate.  As  an  example  we  plotted 
the  variation  of  Ty*  as  a  function  of  wavelength  for 
our  sample  system  (Sec  3.3).  See  Fig.  17  bottom 


Bit  Rate  System  —  Overdense. 


In  this  section  the  optimal  average  system 
throughput  for  the  underdense  burst  is  given.  This 
is  the  maximum  average  number  of  bits  per  unit  time 
for  the  'best'  transmission  rate  R  (and  other  system 
constants)  assuming  only  overdense  bursts  are 
present. 

This  optimal  throughput  is  found  by  first  deriving 
the  dependence  of  the  average  throughput  (Eg. 
3. 4. 3-1)  on  the  bit  rate  R,  and  solving  for  the 
optimal  bit  rate  R*  that  would  yield  the  desired 
optimal  throughput  Tq*  •  The  solution  is  done  with 
respect  to  the  normalized  bit  rate,  RN0,  which  yields 
a  result  independent  of  a  particular  choice  of  system 


values. 
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For  the  given  system  Pbmax,  N0»  R*  a  and  °0  are 
specified.  With  a  reference  to  Eq.  (3. 4. 1.2-3)  and 
the  accompanying  disscussion  we  rewrite  RN0  —  the 
bit  rate  normalized  with  respect  to  Raaxo'  the 
maximum  bit  rate  achievable  for  the  given  system 
using  overdense  bursts  —  and  related  parameters: 


Qt  “  9  uo'  R  NO 


Q  m«o“max{gr,gIO} 


(3. 4. 3. 2-  1  ) 

where,  quo  «  1016  and  q^  =  io14 

a  and  C0  are  defined  in  Eq.  (2. 1.5-2). 

No  is  defined  in  Sec.  2.4 


We  note  from  the  above  that  for  qT  <  g^  (or 
equivalently  RN0  <  .1)  where  qaino  -  qio  and  the 
other  for  q<j  >  qLo  ( or  equivalently  RN0  >  .1)  where 
Smino  "  <3t*  In  addition,  the  value  of  Q  (Eq. 
(3. 4. 1.2-3))  varies  with  the  RN0  as  well  since  Q  is 


a  function  of  q^^  -  see  Eq.  (2. 1.6-2).  For  the 
case  where  %o  >  .1  (q<r  >  qj/j)  we  have 

Q  *  *5 (quo)* 5  RN0  /  (1-rN0) 

For  the  case  where  Rjjo  <  • 1  (<2t  >  Slo)  we  have 

Q  -  .5(qLo) *5  /.9 

Using  the  above  relationships  in  Eq.  (3. 4. 1.2-3)  to 
express  the  average  overdense  burst  duration  as  a 
function  of  RN0  we  get: 

For  Rno  i  0.1  : 

li- ^qLo^u°a'  R^yjjr~  1  "V^T”  1  +  C0S ''(lOR^-cos  ~1{Rno)\ 

(3. 4. 3. 2-2) 

For  Rno  >0.1  : 

h  “  -95Vg  loQuoCL'  Rno{  J  ~rr~  1  -  cos'1  ( *„,<,) 

V  V  Rno 

The  average  number  of  bits  per  overdense  burst 
is  by  definition 

Hbo-R-Tb  (3. 4. 3. 2-3) 


and  can  be  expressed  in  terms  of  the  normalized  bit 
rate  using  Eq.  (3.4. 3.2-1)  for  R  as  a  function  of  RN0 


EJU'JUUV  LIU  «  W  W  W  WW.7  W.  V*  YJ  TV  » ^  « JTTV  V.1  •  JHJ  TV  TX  TX  V>’  ’>  V  WX  V"^r.,r 
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and  Eq.  (3. 4. 3. 2-2)  for  the  average  overdense  burst 
duration. 

Since  the  throughput  Tq  according  to  Eq.  (3. 4. 3-1)  is 
the  product  of  the  arrival  rate  and  the  last  equation 
i.e. 


T  o  “  ^  o  "  ^  «o 


(3. 4. 3. 2-4) 


we  need  to  express  the  arrival  rate  of  overdense 
bursts  Aq  as  a  function  of  rno  .  We  start  by 
employing  Eq.  (2.3-2)  with  q1  set  to  qmino  (the 
lowest  overdense  q  as  determined  from  the  choice  of 
Pmin  bit  rate)  and  q2  set  to  quo  (the  physical 

maximum) .  The  dependence  of  qnino  on  RN0  is  then 
used  (disscussed  above  after  Eq.  (3. 4. 3. 2-1)) 
yielding  the  following  arrival  rate: 


f  .99- 

A‘ 

rno<.\  \ 

A0-{ 

1 

^.01  -  A' 

1 

R-N0>  •  1  J 

(3. 4. 3. 2-5) 

u  1 

_  R  NO  _ 

where , 

M.  3-*sin*(R*£)2 

A  * - 

Q  LO 

quo  *  1016  and  qw  -  1014 


n"  p  i1  g  iu  ^  gnp;  pro  i*  x?  n  v."  g»  iw  v*v^v,v5v»7»v  *  *v^i  «vth*v’.»«»w 
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By  direct  substitutions  from  the  above  equations 
we  express  the  throughput  using  overdense  bursts  as 
a  function  of  normalized  bit  rate  and  system 
parameters  as: 


.33- 10 


1B  PTCTGM{3tsintsec3t-  Xs 


No- g’'{P^a.)- D- Rt 

~l*cos  (iO/?vo)-c os  [p No)^j 
1  Jr-l-cos-*(Ri/o)\ 

K-no  } 


rno<.i 


Rno  R* 

1  —  R„a 


Rmo  ^  •  1 


(3. 4. 3. 2-6) 


where , 


Rt  and  * 


PT 

'  gR 
X 

D 

NO 


9”1 (Pbmax) 


idealized  antenna  beamvidth. 
geometric  parameters  detrmined 
by  link  distance  (from  section  2.2) 
transmitter  power. 

transmitter  and  receiver  antenna  gains, 
carrier  wavelength. 

diffusion  coefficient  of  the  atmosphere, 
power  spectral  density  of  the  noise. 
(Sec.  2.4) 

minimum  of  bit  energy  to  Nq  ratio. 
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Note  that  the  equation  is  written  in  two  parts;  the 
first  incorporates  all  the  parameters  of  the 
particular  system  and  the  second  part  reflects  the 
variation  in  throughput  with  choice  of  bit  rate.  The 
above  throughput  (  multiplied  by  2.55)  for  our  sample 
system  (Sec  3.3)  is  plotted  in  Fig.  18. 

He  would  like  to  find  from  the  above  equation 
the  optimal  normalized  bit  rate  R^o*  which  yields  the 
maximum  throughput  Tq*  .  This  optimal  normalized  bit 
rate  for  any  given  system  is  found  numerically  to  be 
(see  Fig.  18) : 


RN0*  -  .1 

or  equivalently  the  optimal  bit  rate, 

R*  ■  -1  Rmaxo  (3.4. 3. 1-7) 

The  result  is  intuitively  obvious  since  Rjio**1  is  the 
maximum  bit  rate  that  still  includes  all  overdense 
trails  i.e  maximum  arrival  rate.  For  a  higher  bit 
rate  we  have  a  smaller  arrival  rate  that  reduces  the 
throughput.  K  smaller  bit  rate  will  not  change  the 
overdense  arrival  rate  but  will  diminish  the 
throughput . 


i 
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Using  the  optimal  bit  rate  for  overdense  bursts  we 
get  the  maximum  overdense  throughput  To*  by 
substituting  Eg. (3. 4. 3. 2-7)  into  Eg.  (3. 4. 3. 2-6) 
yielding: 


The  above  expression  for  the  maximum  throughput  T0* 
allows  us  to  guantify  the  affects  of  various  system 
parameters  and  constraints  (wavelength,  power,  link 
distance,  maximum  bit  error  rate  etc.)  on  the  best 
throughput  using  a  constant  bit  rate.  Using  our 
sample  system  (Sec  3.3)  with  a  wavelength  of  6  m.  (50 
KHz)  we  get  R*  -  252  kb/s  and  Tq*  ■  8688  b/s. 


In  this  section  we  combine  the  result  of  the 
last  two  sections  to  present  the  total  throughput  and 
optimal  throughput,  since  we  used  differently 
normalized  variables  for  underdense  and  overdense  we 
shall  now  employ  a  single  normalized  variable;  the 
one  corresponding  to  the  overdense  case. 


Recalling  previous  definitions: 
For  underdense  we  used 


** 


R 

R  max 


( 


Q  min 

<7  u 


) 


where 


R 


mmx 


C  uQu 

N  oQ  bntax) 


(3 .4 .3 .3  -  1  ) 


the  maximum  bit  rate  achievable  in  a  given  system 
using  underdense  bursts.  Cu  and  N0  are  defined  in 
Eg.  (2.11-2)  and  Eg.  (2.4-2),  respectively. 


For  overdense  we  used 

R  -  * 

K  NO  R 

maxo 

where 


CoJcT 


<7  vo 


yfe-N0-g-'[Pbnax) 


maxo 


(3. 4. 3. 3-2) 


I 


V| 


$ 


where,  qyo  ■  1016  and  qLO  *  1014 

a  and  Cq  are  defined  in  Eq.  (2. 1.5-2) 
Nq  is  defined  in  Sec.  2.4 

From  the  above  we  can  write 


^-.-A  AAA.  (33.3S9/^MCJ#)i 


s -  6.446 • e 


(3. 4. 3. 3-3) 


Rn  m  s  •  Rt 


(3. 4. 3. 3-4) 


Note  that  R^-l,  the  maximum  underdense  bit  rate, 
corresponds  to  %q  *  1/s . 


Using  Eq.  (3.4.3. 1-4 )  for  the  underdense  throughput, 
with  R)f  expressed  by  the  last  equation,  and  Eq. 
(3. 4. 3. 2-6)  for  the  overdense  throughput  we  get  for 
the  total  throughput  T: 


id  PtGtG.S 20sin0sec2^- Is  ,  , 

T  -  .509  •  1  O'"  ■  r  T  - - (T„  +  T0,*T0,) 

N  0 ' 9  {•  bmax  )'D  Rt 


(3. 4. 3. 3-5) 


where 


T¥m exp 


("Si •  !«(***,)]•  if  (J- a..) 


.,  see +R{ 
to>m  .1401 


p  V **co#*,(IOi**o)"c#,',(Jl»s)]‘£,(-l"Jl*o) 


T02  -  .649  •  sec$ 


&N0  R  wo 


'  1 

— —  1  -  cos*1 
Rno 


(  ^  A/o)  '  U  [ft  NO  ~ 


£  AA,  (33.359/ ls*#eJ<) 

s»  6.446  -ev  -sec$ 


/I  for  x>0\ 
\0  for  x  <0/ 


is  the  unit  step  function. 

All  other  parameters  are  defined  in  previous 
sections . 


For  a  variable  bit  rate  system  the  bit  rate  is 
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changed  so  as  to  maintain  a  reasonably  constant  energy 
per  bit.  Clearly  it  implies  a  constant  BER,  Pj-,,  for 
all  transmitted  bits.  Considering  for  example  the 
underdense  burst,  if  apriori  knowledge  of  the  burst 
were  available,  we  could  divide  the  received  power 
curve  into  equal  energy  portions  and  compute  the 
resulting  individual  bit  durations  before  we  transmit. 
Under  such  conditions  all  bits  will  possess  exactly  the 
same  energy  and  hence  the  same  BER.  The  total  number 
of  bits  transmitted  for  that  burst  will  simply  be 
proportional  to  the  total  area  under  the  received  power 
-versus-  time  curve.  Since  we  don't  have  an  apriori 
knowledge  of  the  burst  but  rather  an  instantaneous 
(within  a  negligible  delay)  indication  of  channel 
strength  we  can  only  evaluate  the  necessary  duration 
for  the  current  bit  and  transmit  it  as  such.  During 
the  duration  of  the  transmitted  bit  the  power  received 
drops.  If  the  bit  duration  is  small  the  resulting  drop 
in  bit  energy  is  negligible  and  our  assumption  of 
constant  bit  energy  is  not  violated.  As  time 
progresses,  however,  the  bit  duration  increases  (power 
received  decays)  and  the  drop  in  bit  energy  relative  to 
the  assumed  constant  bit  energy  becomes  appreciable  and 
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finally  intolerable.  This  fact  limits  the  maximum  bit 
duration  which  in  turn  restricts  the  duration  of 
transmission  or,  equivalently,  channel  duration. 


.5.1.1 


Transmission  Duration  for 


Variable  Bit  Rate  —  Underdense 


We  start  by  imposing  our  fidelity  constant,  k, 
for  the  normalized  bit  energy  —  the  ratio  of  the 
real  bit  energy  to  the  desired  one.  For  our 
derivation  we  will  assume  k  >  .9  i.e.  no  more  than 
10%  drop  in  bit  energy  shall  be  tolerated.  We 
further  define: 


Pb 

Tbi 

R(t) 

To 

N 

T(m) 


Assumed  constant  BER  for  all  transmitted 
bits. 

Duration  of  the  i^h  bit 
Instantaneous  bit  rate 
Duration  of  the  first  bit,  1/Rq 
Last  bit 

Total  duration  of  the  first  m  bits; 
also  starting  time  of  the  (m+l)th  bit: 


i*  i 


! 

I 


wowWKWOBTWOTOWWWPTOPTiPTionpnonpnBnDnDrgygwBimroniTCroiOTP  gCTCTenrosF^wBCTwrarTqngron 


I 

I 
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Tt 


/?(0“ 


T(«-l) 


Duration  of  transmission  (burst  duration) 
T(N) 

Max  (Tjjjj/B)  for  a  given  k. 


In  order  to  maintain  a  constant  BER,  P^,  we  have  for 
the  variable  bit  rate 

cvq2 


-e*"'-  Rn-e',/t 


No  g"{P>) 

and  the  n^*1  bit  duration  is  given  by 


(3.5.1  .1  -  1) 


T^-ToS  ' 

where 

T0  -  1/R0  -  N0  g“l  (Pb)/Cu  q2 

For  k  >  .9  we  have  for  the  last  bit,  N: 

Tb„<[\  .5->/2. 25-6(1  -  It)]  -  B-c  B 

Combining  with  eq  (3. 5. 1.1-2)  we  have 


(3.5.1  .1-2) 


(3.5.1  .1  -3) 


z±»_z_ 

B  ’  B 


0  e  *  <c 


(3.5.1  .1-4) 


Since  the  transmission  duration  is  given  by 

Tt  "  T(N-l)  +  TfcN  we  have  from  the  last 

3  eqns. 


7\<B-[c  +  ln|  C-A  I] 


(3.5.1  .1-5) 


for  k  ■  .9,  c-  .2  we  have  from  the  above  T^/B  -  .2  + 
In  (R0B/5) . 

For  all  practical  systems  ve  can  safely  use 


7,  -  £  •  In 


R0B 


(3.5.1  .1-6) 


where  Rq  Is  defined  In  Eg.  (3. 5. 1.1-1)- 


Using  the  results  from  the  previous  section  on 
transmission  duration  we  can  find  the  approximate 
number  of  transmitted  bits  by  integrating  R(t)  from 
time  0  to  time  Tt: 

NBU“R0B"cUB<32/9_1(pb)N0  (3. 5. 2-1) 

where  R0  is  the  variable  bit  rate  at  time  0. 

and  since  B  and  q  are  independent  the  average  of  Nyy 

over  the  joint  space  of  q  and  B  yields: 


UV  -  \ 


C  ulfq1 
g"{Pt)-N< 


(3. 5. 2-2) 


The  derivation  of  the  above  equations  is  in  section 


We  shall  assume  that  the  bit  rate  shall  vary 
between  RminV  and  RmaxV*  Bearing  in  mind  that  for  a 
given  bit  rate  R  we  have  a  corresponding  q  such  that 


C  uQ  “N0'g  ( P  b)’R  Q  l<  Q  <  Q  u 


As  a  result  we  have 


(3. 5. 2-3) 


Q  minV 


Q  maxV 


'  No-g~l[Pb) 


'Wo'g-’tM 


(3. 5. 2-4) 


(3. 5. 2-5) 


RminV  must  be  chosen  such  that  qminv  >  Sl*  the  physical 
minimum. 

RmaxV  oust  be  chosen  such  that  qmaxV  K  9u»  physical 
maximum. 

Note  that  because  our  empirical  arrival  rate  relation¬ 
ship  (Sec.  2.3)  is  valid  for  q>1013  we  shall  set  RminV 
such  that  qBinV  >  1013* 

To  get  the  most  throughput  out  of  our  system  we  set 


'*  "V  •  »*  *>  V** 


WWW 
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RmaxV  to  its  maximum  value  i.e.  corresponding  to  q=qu# 
We  thus  have  set  qmaxV^u  and  using  it  in  the  last  two 
equations  we  have 


'‘♦1 


ll_ 


Q  minV  Q  minV  I  R  minV 
Q  maxV  Qu  V  RmaxV 


V  Rnv 


(3. 5. 2-6) 


where  RjnaxV  i®  given  by  Eq.(3.5.2-5)  as 

_  CuQu 

u0g  '{p b) 


(3.S.2-7) 


where,  Cg  is  from  Eq.  (2.1. 1-2) ,  Nq  is  from  (2.4-2) 
resulting  in  Eq. (3.5.2-13) . 

Note  that  RmaxV  “  Rmax  the  maximum  underdense  for 
constant  bit  rate  system  if  we  set  Pb^Pbmax*  a  direct 
result  from  Eq.  (3. 4. 3. 1-1). 


Evaluating  E{q2)  for  Eq.  (3. 5. 2-2)  in  the  range  of 
qminV*  qu  i®  done  using  eq  (2.12-2) 


— L>qp-i  g 3-p (  ~y3  P) 

H  (3-  Pym'nvqu  (  J  _  y  P'  1  j 


where  p«1.6  and  (with  Eq.  (3. 5. 2-6) 


Q  minV  r~Z 

y -——m\IRnv 
Qu 


Rewriting  we  get  for  the  last  equation 


(3. 5. 2. -8) 


I 


Si 


I 


1 


3 

h 


i 


3 


.V 
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— 2  *6  2!^^  ^  Ar^ ) 

<7  “7— T<7(/ 


1-4""  (l-K*3,) 


(3. 5. 2. -9) 


Substituting  into  Eq.(3.5.2-2)  we  get 


N 


C  u~E  .6 


yv'  .  1 

g 


[P b)-N0  1  -4Qt/  (1  -Rjv) 


(3.5.2-  10) 


For  throughput  expression  we  need  to  have  the 
arrival,  Ay,  expressed  as  a  function  of  our  variable 
Rjjy.  From  Eg.  (2.3-2)  we  have 


av-^-Z-  (i-2i=r)-i2f 


g  mmv  \  g  u  )  g  u  \  Jr 

where 

A'v  -  19. 1667-  <j>sin<p(RRZ)2 


(3.5.2-11) 


NV 


The  expression  for  the  throughput,  Tyv  is  dictated 
by  Eq.  (3. 4. 3-1)  and  is  arrived  at  by  the  product  of 
Eq.  (3.5.2-11)  and  Eq.  (3.5.2-10)  above.  Using  the 
definitions  of  Eq.  (2. 1.1-2)  and  Eq.  (3. 5. 2-7)  we  get 


T  vv  =  6.54821  •  10 


-20 


p  T  0  T  C  *  £  2  0  s  i  n  ( <p  )  s  e  c 2  ( <p  )  e  x  p  {  -  V : 

\  A  s«c  4>  / 


RTKog~'(Pt)D 


A 


a 


1 


■vrvcviva 


8 


I 


! 


i 


s 


HE 

Vj» 

Sv 

2 


& 


.>/* 


:-  1 


1  .  u  ^  NV  R  NV , 


NV 


(1  -Ri 


3  ' 
NV  ) 


100 


(3.5.2-  12) 


where,  (Eq.  (3. 5. 2-7)  ,  Eq.  (2. 1.1-2)  and  (2.4-2)) 


RmaxV  -  2 .5 1 8  •  10 


16 


_  _  „  ,  3  (  33.359  ) 

P  tGtGr  A.  exp  I  —  j 

\.  1  sec  ♦  / 


rnv 

(3. 5. 2-7) 


(.063A.2'3+l),g"'(P  b)'/?r 

RminV/RmaxV  *  RmaxV  defined 


(3.5.2-  13) 

in  Eq. 


RminV  is  arbitrary  yet  if  resulting  Rjjv  is 
small  system  complexity  and  cost  increases. 
K  idealized  antenna  beamwidth. 

RT  and  0  geometric  parameters  detrained 

by  link  distance  (from  section  2.2) 

PT  transmitter  power . 

Giji,Gr  transmitter  and  receiver  antenna  gains. 

a.  carrier  wavelength. 

D  diffusion  coefficient  of  the  atmosphere. 

Nq  power  spectral  density  of  the  noise. 

(Sec.  2.4) 

g”1(Pbmax)  minimum  of  bit  energy  to  Nq  ratio. 


For  our  sample  system  (Sec.  3.3)  we  chose  Rminv=100° 
b/s  i.e  qminV  *  1,1  1°13'  g_1(Pfc)=9  we  *?et  RNV=* 011558 
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and  Tyv  about  1150  b/s. 


If  we  assume  the  constant  bit  energy  of  the  variable 
rate  system  tobe  the  same  as  the  minimum  bit  energy  of 
the  constant  rate  system  then  for  our  sample  system  ( 
underdense  only)  we  improved  the  throughput  from  214 
(Sec.  3. 4. 3.1)  to  1150  -  a  ratio  of  5.3. 


3.5.3  Throughput  for  Variable  Bit  Rate  Svste 


The  maximum  power  point  of  an  overdense  burst  PB 
is  well  aproximated  by  (Eq.  (2. 1.5-2)  neglecting  k) 


P  max  Co 


(3.5.3-  1  ) 


We  recall  that  for  our  variable  rate  system  the  rate 
varies  as  the  power  scaled  by  the  bit  energy  Nog_1(Pb) 
This  variation  in  time  can  be  approximated  as  a 
triangle  whose  base  lies  on  time  axis  from  the  origin 
to  the  point  given  by  aq  and  whose  top  vertex  is  at 
Pmax/Nog_1(Pb)- 

Since  the  number  of  bits  sustained  by  this  overdense 


i 


burst,  Ngv»  is  theoretically  equal  to  the  area 


underneath  the  rate  tine  function  we  have 


>,  1  Pmmx 

N  iv  —7—-. '  a- q- 


,  /-  _  3/2 

1  c0a  3/2 


2  N0g‘l{Pb)  H  2{eN0g'l(Pb) 


(3. 5. 3-2) 


Averaging  over  the  ensamble  of  overdense  bursts  using 
Eq.  (2. 1.6-2)  with  qnino“<JLO“1014  (9L0/9U0  "  *01) 


(3. 5. 3-3) 


which  yields  for  average  number  of  bits 


^  nv  “ 


1  C0a3/3 

4^A/og-‘(/>b)' 


(3. 5. 3-4) 


and  the  arrival  rate  is  given  by  Eq.  (2.3-2) 


3-*sin  <7 


f1-— ) 

V  Quo) 


(3. 5. 3-5) 


Combining  all  of  the  above  for  the  throughput 
expression  (varible  rate  and  overdense  bursts)  we  get 
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ov 


94.686- 10 


-JO 


/‘rCrC^fsinfsec3#-  Is 


(3. 5. 3-7) 


No-g"[rb)-D-RT 

For  our  sample  system  (Sec.  3.3)  we  get  a  throughput  of 
26.7  kb/s. 


If  we  assume  the  constant  bit  energy  of  the  variable 
rate  system  to  be  the  same  as  the  minimum  bit  energy  of 
the  constant  rate  system  then  for  our  sample  system 
(Sec  3.3)  (assuming  overdense  only)  we  improved  the 
throughput  from  8688  b/s  (Sec.  3. 4. 3. 2)  to  26.7  kb/s  - 
an  improvement  ratio  of  3. 


3.5.4  Comparison  of  Throughput  for  Constant  and 
Variable  Bit  Rate  Systems 


In  this  section  we  shall  apply  the  derived 
relationships  in  the  previous  sections  to  a  practical 
system.  The  various  parameters  of  interest  are 
computed.  The  throughput  of  the  system  for  the  optimal 
constant  bit  rate  is  computed  and  compared  with  that  of 
throughput  of  the  system  using  variable  bit  rate.  This 
comparison  will  demonstrate  the  improvement  in 
throughput  for  variable  bit  rate.  We  start  by  copying 


the  parameters  of  the  system  from  Sec.  3.3: 


L 


sec(^) 

At-** 

PT 

G<ji 

gR 


■  1000  km  link  distance.  From  section  2.2 
we  have  for  the  next  two  parameters: 
*18.5  and 
*514  km. 

-  1000  w 

-  10  dB. 

*  10  dB.  And  for  antenna  beamwidth  we  use 

=  45.84°  (.8  rad) 


A.  =  6  m.  (50  MHz) 

D  *  8  m2/s 

Pbmax  “  10"5  Max.  BER  for  constant  bit 

rate  system. 

Using  BPSK  modulation  we  get: 
g’1(Pbmax)  -  9 


W/Hz 


See  discussion  above  Eq.  (2.4-2). 

N0  -  4.89  10“20  at  50  MHz . 

For  comparison  the  BER  for  variable  rate  system,  P^, 
is  set  equal  to  Pbmax 

RminV  “  1000  b/s  (minimum  variable  rate  allowed) 
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The  improvement  we  get  for  underdense  burets  is  arrived 
at  under  the  assumption  that  the  optimal  rate  is  used 
for  the  constant  bit  rate  system  and  for  the  variable 
bit  rate  system  some  Rninv  is  assumed.  In  addition,  we 
assume  that  the  Pb*Pbmax  i.e.  the  constant  bit  energy 
(variable  rate)  is  assumed  equal  to  the  minimum  allowed 
bit  energy  (constant  rate  system) .  We  thus  divide  the 
relevant  equations:  Eq.  (3.5.2-12)  and  Eq.  (3. 4. 3. 1-6) 
resulting  in  the  improvememt  lw 


"  J  (i  -*&) 

where, 


RmaxV  -  2.518  •  10 


(.06312  3+ 


(3.5.2-  13) 


The  above  is  evaluated  using  the  values  of  our  sample 
system  where 

RNV“Rminv/Rmax  ™  1000/86500  ■  .01156 
resulting  in 


Itt  ■  5.38  or  about  7  dB 


The  relevant  equations  to  be  used  for  comparison  of 
throughputs  are  Eq.(3.5.3-7)  and  Eq.  (3. 4. 3. 2-8) 
resulting  in  an  improvement  of 


I0  “  T0v/T*o  "  3 


Disscussion 


|  We  shall  compare  now  the  improvement  in  throughput  of 

adaptive  rate  relative  to  fixed  rate  systems  for 
|  underdense  in  terms  of  its  components;  the  arrival  rate 

and  the  average  number  of  bits  transmitted  per  burst. 

*  For  fixed  rate  (underdense  bursts)  we  have  the  folowing 

results : 

|  Fixed  Rate: 

Umax  m  ® ® ®  Kb/ s 

\  Since  the  normalized  optimal  bit  rate  is».07  we  have 

'  R*“  •  07  Rgux  m  8  •  05  kb/ S 


For  our  system  we  also  have 


"i, 


»* 


1 


W 


■o\ 

M 


U.1 

M, 


i  i  ■'«  i  •  ft  (  1^  ftj*Ll « S  |*|  ft«>  «*l  a’i  »**  ■•> 


5  -  .264  sec 


Cv- 3.8-  10' 


From  eq  2.4-2  we  have: 


N0  -  4.89  10“20 


For  this  rate  we  can  compute  some  parameters  for 


general  interest: 


the  average  number  of  bits  per  underdense  burst  is 
(from  Eq.  (3. 4. 3. 1-2))  1,846  kb/s  per  burst 


<3min  "  .26457  1014  e/m  (  from  eq.  3. 4. 1.1-1  set  R  to 
R*) 


Which  yields  an  arrival  rate  of 


•  11$  bursts  per  sec  or  equivalently  an  average 


inter-arrival  time<-8.6  s. 


The  product  of  the  last  two  results  yield  a  throughput 


of  214  b/s 


variable  Rate 


In  general  the  ratio  of  the  bursts'  arrival  rates  of 
variable  bit  rate  to  fixed  bit  rate  system  is  given 
by 

Ay/A*  -  .36  {  (Rmax/Rminv)  * 5  "1  >  (Eq.  (3. 4. 3. 1-3) 

and  Eq. (3.5.2-11) 

which  in  our  example  yields  2.9. 

The  average  number  of  bits  per  burst  using  variable 
rate  is  from  Eq.  (3.5.2-10) 

3.328  kb/s  per  burst . 

The  ratio  of  improvement  for  the  average  number  of  bits 
per  burst  is 
3.328/1.846  «  1.8 

whereas  the  improvement  in  arrival  rate  is 
2.9  yielding  a  total  improvement  of  5.3- 
It  is  interesting  to  note  that  the  larger  improvement 
is  due  to  the  arrival  rate.  The  arbitrary  set  RminV, 
the  minimum  rate  allowed  in  adaptive  rate  system,  sets 
a  lower  threshold  of  electron  line  density,  qminv» 
whic  utilizes  bursts  that  are  precluded  for  the  fixed 
rate  system  with  its  higher  electron  line  density 
threshold  q^in* 
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4  DERIVATIONS 


This  chapter  supplements  the  derivations  outlined  in 
previous  chapter  for  the  results  that  have  thus  far  been 
presented. 


With  reference  to  fig.  6  and  the  definitions  in 


section  2.2  we  have  the  following  relationships: 
P-L/4R, 


y-rp 


x  -  2/?, sin 


(4.1-1) 


(4.1-2) 


(4.1-3) 


y  =  xsin(i&)  =  2/?,sin2(t>)  (4.1-4) 

Since  it's  always  true  (L  <  2000  km)  that  (L/4Re)2  «i  we 


have 


sin2(*>)-*>2-  -jr- 


(4.1-5) 


By  the  Law  of  Cosine  we  have: 

h2  +  x2-2xhcos(y)  = = /?2  (4.1-6) 

We  substitute  for  x  and  y  and  smJ(*o  from  eqns.  (4.1-3), 

(4.1-2),  (4.1-5)  respectively  in  eqn.  (4.1-6)  and  arrive 

at  the  stated  result  in  section  2.2  after  some  trivial 


S 


algebra  : 


^T-^«“V(/l+l2/8/?.)2  +  i2/4 


(2.2-1) 


From  the  figure  we  also  have: 


sec2($) 


Rj 

(h  +  y)7 


which  when  combined  with  eqs. (4.1-4),  (4.1-5)  and 

(2.2-1) 

yields: 


sec2(0) 


1  ♦ 


(20^)J 


(2.2-2) 


and  thus  completing  the  derivation. 


4.2  Derivation  of  Average  Burst  Duration. 
Constant  Bit  Rate 


This  derivation  is  for  the  underdense  case.  We  start 
with  equation  (3. 4. 1.1-1)  as  the  expression  for  tB> 
Since  we  want  to  average  over  the  joint  space  of  q  and  B 
which  are  statistically  independent.  The  joint  pdf  is 
given  by  the  product  of  the  individual  pdfs.  Because  B 
acts  as  a  coefficient  in  the  expression  for  tB  the  two 
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integrals  are  separable.  (The  individual  pdfs  for  q  and 
B  are  from  equations  (2. 1.2-2)  and  (2. 1.3-2, 3)  respec¬ 
tively. 

tKB-2,cB*-Q  C'  {iniq/q^Y-q-'dq  (4.2-1) 

•'i.  H 


The  above  integral  can  be  solved  by  employing  the 

following  substitutions: 

z  -  (p-1)  In  (q/qmin) 

q  *  <3min  e*/(P"l) 

dq  *  tqmin  /(P“D]  e*/<P-U  dz 

lnn(q/qn,in)-C2/(p-l)  ]n 

Reduces  the  above  to  a  single  integral  of  the  form 


2  nB* 
(P-D" 


1 

l-x 


ln(  I  / x) 


zne 


Zdz 


(4.2-2) 


where 


(  QmiuY'1 

~N0-R’g-'{Pbnaxy 

\  Qu  J 

CuQu 

which  can  be  solved  recursively  by  parts,  yielding 


2  "IP  n1 

(p-1)*  l-x 


i-0 


ln‘(x)\ 


Setting  n«l  and  n«2  in  the  above  yields  the  necessary 
equations  for  the  expected  value  and  variance  of  tg 
although  some  wading  through  algebra  is  necessary  for  the 


Since  the  roots  of  the  overdense  power  received 
equation  cannot  be  solved  for  directly  we  use  a  second 
order  approximation.  In  fact  each  side  of  the  overdense 
'mountain*  is  fit  with  a  parabola.  The  parabolas  are 
expanded  about  the  maximum  power  point.  For  the  left 
parabola  (which  approximate  the  left  side)  we  use  the 
smaller  root  and  subtract  it  from  the  larger  root  of  the 
right  parabola  (which  approximates  the  right  side  of  the 
power  equation) .  This  difference  is  the  burst  duration 
for  the  overdense  trail.  We  start  by  copying  eq. 
(2. 1.5-2) 

P(t)-CoyJ(k  +  t  <aq-k  (2. 1.5-2) 


where. 


V 

$ 
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C0  -3.16628610  3 fD ——■  -  ■ 

Rt 


k  -  .  1  05625/Z) 


a  -  7. 143143 10 
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We  need  to  solve  for  the  roots  t2,  t^  of  the  above  eq. 
when  P  is  set  to  Pmin  since  their  difference  is  our 
transmission  interval  -  burst  duration.  To  do  so  we 
first  modify  the  above  eq.  as: 


r,lnlr. 


k/ aq  <  r  <  1 


(4.3-2) 


or  equivalently 

rln(r)+  P7N  =  0 

where 

Pn  -  P  /  Co  (aq) 1/2 
r*(k+t)/aq 

The  above  will  yield  the  roots  r2,  for  P=Pmin‘ 
Note  that  the  burst  duration  is  defined  as 


(4.3-4) 


tB*st2-ti=aq(r2-r1) 


(4.3-4) 


We  are  now  ready  to  approximate  r2-ri  where  r^  is  the 
•left'  root  and  r^  the  'right'  root.  Note  that  the  above 
eq.  has  a  maximum  P^  of  l/(e)l/2  at  r=l/e. 


To  approximate  we  expand  ln(r)  about  r-l/e  and 

substitute  the  first  two  terms  for  ln(r)  in  the  last  eq 

resulting  in:  r2  -  (2/e)r  +  (l/e)Pji2  -  0 

Solving  and  using  the  smaller  root  we  thus  approximate  r^ 

as: 


To  approximate  T2  we  expand  rln(r)  about  r=l/e  and 
substitute  the  first  three  terms  in  eq  (4.3-3)  and  solve 
for  the  larger  root. 

(4.3-6) 

(4.3-7) 

A  comparison  of  the  above  with  an  iterative  solution  of 
the  real  (without  approximation)  roots  of  the  eq  (4.3-3) 
and  their  difference  allows  us  to  replace  the  .88 
constant  in  the  last  eq.  with  .95 

r2”rl  “  *95  root  [1-e  P  sub  N  super  2] 


i 


Since  the  power  for  which  we  want  to  find  the  roots  is 
Pmin  fr°®  eqs.  (4.3-3)  and  (4.3-4)  we  have  the  desired 


* 


! 


1 


approximation  for  tg 


n 


i 

5 

I 

g 

I 

i 
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s 

4m 

IS 


t,m  .95a-  <7-  a/  1  - 


_ u  e  1 

C*a  q 


(4.3-8) 


which  is  the  result  stated  in  (3. 4. 1.2-1). 

Note  that  q>qmino  ■  max  (St*  SLo)*  If  3t  (which  is  set 
by  Pmin  and  system  parameters  reflected  in  Co  and  a  is 
less  than  q^o  than  the  minimum  overdense  electron  line 
density,  q^ino*  is  set  to  the  Physical  minimum  overdense 
electron  line  density  qj^.  On  the  other  hand  if  qT  is 
bigger  than  q^  than  the  overdense  trails  with  electron 
line  density  between  qj^,  and  q^  have  zero  burst  duration 
and  are  not  "seen"  by  our  system. 


For  most  practical  systems  which  utilize  underdense 
bursts  P]nin  is  such  that: 

q?  «  qLo  ^mino  *(lLo  and  is  approximated  as  aq, 

(q>qLo  )  and 

E(tB)  -.95a  E(q)  -  a(quo  Qlo)  V2  -  a  10^. 

Derivation  of  CDF  and  PDF: 


r 
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From  eq.  (3. 4. 1.2-1)  we  note  that  tg  is  monotonically 
increasing  with  q  such  that  for  any  tg<t^  we  have 
correspondingly  a  q<qj..  For  tg  ■  t^  we  have  from  eq. 


(3.4. 1.2-1)  : 


(4.3-8) 


where  r  =  agr/2  =  (e/2)(/>min/C0)2 

From  the  above  discussion  there  is  a  minimum  tj,  timin 
corresponding  to  qmino  and  a  maximum  t^,  tljnax 
corresponding  to  quo>  the  max  overdense  electron  line 
density.  We  can  now  relate  the  following  probability 
events : 


Prob{ta<f1}-Prob{gmino<q<g1} 


(4.3.9) 


We  get  the  CDF  of  q  using  eq  (2. 1.2-2)  such  that  the 
right  hand  side  of  the  last  eq.  becomes 

Q'  tl“(qmino  /  qi)p_1] 

Eq  (4.3-9)  then  becomes 


«i-.»  (4.3-10) 


Substituting  eq  (4.3-8)  into  the  above  eq.  we  get  the 
desired  eq.  (3. 4. 1.2  -6). 

To  get  eq.  (3.4. 1.2-7)  we  simply  take  the  derivative  of 
eq.  (3. 4. 1.2-7) 


•j. 


..  .  ^-icimai  Average  ^nrouanout  ror 

constant  bit  rate  -  Underdense 


Our  first  objective  is  to  derive  eq.  (3. 4. 3. 1-2)  from 
section  3.4.2  we  know  the  number  of  bits  per  given 
underdense  burst  for  a  given  R  and  a  set  of  specified 
system  parameters  is  given  by 

Nru  “  (4.4-1) 

Averaging  Nbu  over  the  ensemble  of  underdense  bursts  (q 
and  B  independent)  is  given  by 

Kbu-R-Tb  (4.4-2) 

rewriting  the  result  in  Eq.  (3.4.1. 1.-3)  which  is  derived 
in  Sec.  4.2 


25 


p  -  1  1  -  x 


[  1  -  x(  1  -  In  x)] 


(3.4.1  .1  -3) 


where 


fq^Y" 


P~ 


5-3.16610~j-  —  ~r-  qu-  10M  p-  1.6 

d  _  ^  ^  mln  \ 

P - 

•»  max  VI/ 

C  vql 

"““iVog-H^n.) 

we  substitue  for  x  and  R  in  terms  of  their  RN 
definitions  and  write 

(RN)P'  =  x  P'  =  (p-l)/2  (p  -1.6) 

R  “  rN  Rmax 

Substituting  these  into  eq.  (4.4-2)  yields  the  desired 
eq.  (3. 4. 3. 1-2) 

nT5  j 

^  BU  ”  1  n  ( ^  )  ]  (3. 4. 3. 1-2) 

P  1  1  -  A  N 

This  eq.  specifies  for  a  given  choice  of  R  the  average 
number  of  bits/burst.  The  rest  is  in  Sec.  3. 4. 3.1 


Rate 


With  reference  to  section  3. 5. 1.1  ve  define  again  the 
following  terms: 


R(t) 


T(m) 


Assumed  constant  BER  for  all  transmitted 
bits. 

Duration  of  the  1th  bit 
Instantaneous  bit  rate 
Duration  of  the  first  bit,  1/R0 
Last  bit 

Total  duration  of  the  first  m  bits; 
also  starting  time  of  the  (m+l)th  bit: 


8 


Duration  of  transmission  (burst  duration) 
T(N) 

Max  (Tfcjf/B)  for  a  given  k. 


We  wish  to  maintain  the  actual  bit  energy  above  a 
certain  fraction  say  90%  of  the  assumed  bit  energy.  This 
will  insure  a  relativly  contant  BER  corresponding  to  the 
assumed  bit  energy.  Since  the  last  bit,  the  Nth  bit,  has 


the  longest  duration  (power  is  monotonically  decreasing) 
the  loss  in  bit  energy  relative  to  the  assumed  level  is 
the  largest,  n  equation  fora: 


Eb  actual  >  *  Eb  assumed 

i.e. 

I  P(T(N-  \))-e~,/,dt>k- P(T(N-  \))-TbN 
J  o 

which  implies 

1  -  kT'bN  >  eJfP(-T'bN)  where  T'bN  -  TbN/B 
and  assuming  T'bN  <  1  we  use  the  first  four  terms  in  the 
series  expansion  of  the  exponential  (about  0)  and  get  the 
root  of  T’bu  yielding 

The  quadratic  equation  in  T'bj|  yields 

.5-72.25-6(1  -  Ac)  -  c 

From  the  above  definitions  of  the  bit  duration  we 

have 

TbN/B  -  tl/Ro  B)  exp  (T(N-1)/B)  <  C 
yet  by  definition  T(N-l)  *  Tb  -  T^jj 
we  have  then  Tb  <  B  [c  +  In  (cR0B) ) 
which  is  the  result  stated  in  eq.  (3. 5. 1.1-5) 


We  now  evaluate  the  number  of  bits  for  a  given  burst 
using  variable  bit  rate  as 


Tiubt 


N 


au 


-  J  R(t)dt- R0bI  1  -exp 


B 


Using  the  maximum  from  eq.  (3. 5. 1.1-4)  i.e.  letting  NB 
«  maximum  Ng  we  have 

nBU  “  R0®  +  1  “1/c 
For  K«.9  c= . 2 

For  all  practical  system  RqB»4 
such  that 


Nbu  -  RoB  -  cu  B  q2/q“1(Pb)  N0 
as  stated  in  section  3.5.2  Eq.  (3. 5. 2-1) 


The  average  throughput  is  proportional  to  the  ensemble 
average  of  NBU.  Since  B  and  q  are  statistically 
independent  we  have 


TT7u 


C  u 


g~'[Pb)-  N( 


which  is  Eq.  (3. 5. 2. -2) 
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5  CONCLUSION  AND  IDENTIFICATION  OF  FUTURE  WORK 


In  this  work  a  complex  stochastic  model  for  the  meteor 
burst  channel  was  used  to  serve  as  the  bedrock  for  the 
analysis  of  two  different  communication  systems.  The  first 
system  assumed  we  are  transmitting  with  a  constant  bit  rate 
while  the  BER  is  time  varying  yet  constrained  to  be  below 
some  maximum  level.  The  second  communications  system 
employs  a  variable  bit  rate  whose  adaptive  behavior  is 
controlled  so  as  to  maintain  constant  BER.  Both  systems 
operate  under  the  control  of  a  communication  feedback 
protocol  that  provides  knowledge  of  channel  arrival  (the 
instant  at  which  the  channel  first  becomes  available)  and 
monitor  the  strength  (power  received)  of  the  channel.  The 
research  herein  focused  on  throughput  and  related 
parameters.  The  average  throughput  was  shown  to  depend  upon 
the  arrival  rate  of  the  bursts  and  the  expected  number  of 
bits  per  burst  as  averaged  over  the  ensemble  of  (random) 
time  functions  that  describe  its  behavior.  In  deriving  the 
expressions  for  the  above  it  was  necessary  to  derive  the 
relationship  between  the  number  of  bits  per  given  burst  and 
bit  rate.  This  in  turn  requires  understanding  of  channel 
duration  -  the  available  transmission  time  -  for  a  given 


channel  (burst) . 
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For  constant  bit  rate  system  the  channel  duration 
decreases  with  increasing  bit  rate.  Since  the  number  of 
bits  transmitted  on  the  burst  is  given  by  the  product  of 
burst  duration  and  bit  rate  a  clear  tradeoff  exists  between 
bit  rate  and  transmission  time.  This,  naturally,  implies 
the  existence  of  an  optimal  bit  rate  that  would  maximize  the 
number  of  transmitted  bits  for  the  burst.  Having  chosen  an 
optimal  bit  rate  for  the  burst  at  hand  we  must  repeat  the 
process  for  the  next  burst  whose  behavior  is  different.  In 
fact,  to  optimize  the  throughput  we  must  change  the  bit  rate 
from  one  burst  to  the  next.  Having  to  change  the  bit  rate 
on  a  per  burst  basis,  however,  requires  knowledge  of  the 
burst  time  behavior.  The  problem  with  this  last  requirement 
is  twofold;  first,  a  priori  knowledge  of  the  burst  is  not 
available  and  second,  estimation  of  the  burst  based  on 
assumed  ideal  behavior  is  unreliable  since  the  bursts  often 
deviate  from  such  assumed  bahavior.  To  counter  this 
limitation  we  derived  the  relationship  between  the  constant 
bit  rate  and  the  expected  average  number  of  bits  per  burst 
since  the  latter  is  directly  proportional  to  throughput.  We 
then  find  the  constant  bit  rate  (for  the  system)  which  will 
maximize  the  throughput.  This  bit  rate  is  fixed  over  all 
bursts  and  therefore  it  is  a  system-optimal  bit  rate. 
Having  found  the  maximum  average  number  of  bits  per  burst  we 
multiply  by  the  arrival  rate  of  burst  and  get  the  maximum 
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throughput  expression  in  terms  of  system  parameters  and 
constraints.  The  significance  of  such  expression  is  its 
utility  in  analyzing/ improving  existing  systems  and  in 
designing  future  constant  bit  rate  meteor  burst  systems. 

In  variable  bit  rate  systems  we  change  the  bit  duration 
(on  a  bit  by  bit  basis)  so  the  bit  energy  remains  constant 
corresponding  to  our  desired  fixed  BER.  The  bit  energy  is 
given  by  the  time  integral  of  the  power  received  function 
over  the  duration  of  the  bit.  Finding  the  individual  bit 
durations  is  tantamount  to  dividing  the  power  received  curve 
into  equal  strips  each  of  area  equal  to  the  bit  energy.  As 
discussed  earlier  a  priori  knowledge  of  the  burst  time 
function  is  not  available  and  estimation  is  at  best 
problematic.  This  precludes  exact  solution  of  the  bit 
duration.  The  best  we  can  do  is  determine  the  bit  duration 
based  on  the  bit  power  at  the  beginning  time  of  the  bit. 
This  may  yield  an  optimistic  (too  short)  bit  duration  if  the 
power  decays  during  the  duration  of  the  bit.  The  energy  per 
bit  for  such  a  case  would  be  below  the  desired  fixed  one. 
If  the  power  is  monotomically  decreasing  there  will  be  a  bit 
whose  bit  energy  relative  to  the  assumed  bit  energy  will  be 
unacceptable.  The  restriction  on  maximum  bit  duration  in 
turn  constrains  theoretically  yet  not  practically,  as  has 
been  shown,  the  total  time  of  transmission  for  the  burst. 
The  number  of  bits  transmitted  during  a  burst  can  be 


approximated  by  the  area  under  the  power  received  for  the 
interval  given  by  transmission  time.  The  throughput  is 
found  by  multiplying  the  arrival  rate  by  the  average  number 
of  bits  per  burst. 

As  it  turns  out  the  throughput  for  variable  rate  system 
is  about  six  dB  higher  for  a  conservative  sample  system  . 
The  interesting  point  however  is  the  fact  that  the 
improvement  of  throughput  for  variable  over  constant  bit 
rate  system  (using  underdense  bursts  only)  was  in  large  part 
due  to  the  arrival  rate.  For  constant  bit  rate  system  the 
high  optimal  bit  rate  precluded  many  bursts  (those  with 
electron  line  density  below  the  minimum  line  density 
threshold)  from  being  useful  resulting  in  a  low  arrival 
rate.  For  variable  rate  systems  since  the  bit  rate  can  be 
quite  low  (though  limited)  implies  a  lower  level  of  cutoff 
for  useful  bursts  and  higher  arrival  rate.  The  tradeoff 
lies  in  the  fact  that  a  low  threshold  rate  for  the  adaptive 
rate  increase  system  complexity. 

The  improvement  in  throughput  however  must  be  kept  in 
perspective.  Our  analysis  assumed  that  the  constant  BER 
under  variable  rate  system  is  same  as  the  worst  BER  for  the 
constant  bit  rate  system.  Consequently,  given  N  transmitted 
bits  the  probability  of  all  being  error  free  is  lower  for 
the  constant  bit  rate  then  for  the  variable  bit  rate  system. 
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The  higher  throughput  for  overdense  relative  to 
underdense  bursts  Is  a  natural  outcome  of  Its  longer  average 
burst  duration.  The  tradeoff,  however,  is  in  the  much 
longer  inter-arrival  time  of  the  overdense  bursts  compared 
with  underdense;  a  ratio  of  15:1.  We  see  then  that  for 
short  urgent  messages  where  waiting  time  is  crucial  if  we 
are  given  a  constant  bit  rate  system  it  is  advantageous  to 
operate  the  system  at  the  optimal  rate  corresponding  to 
underdense  burst  utilization. 

In  this  analysis  many  parameters  considerations  and 
variables  interplay.  Amongst  them  complexity  of  model, 
empirical  values,  transmitter  power  antenna  design  link 
geometry  and  geographic  location,  frequency,  noise  behavior 
modulation  bit  rate  and  BER.  All  of  the  above  have  been 
integrated  and  the  derived  closest  form  relationships  are 
stated  in  the  most  general  form  incorporating  all  possible 
variables  to  facilitate  their  usage  as  design  and  analysis 
tools.  Furthermore,  these  expressions  define  and  better 
delineate  the  tasks  for  future  analysis.  Some  of  the  issues 
that  need  be  investigated  are  coding,  packet  design, 
modulation  protocol  development  and  networking.  Each  of 
these  requires  extensive  research  and,  clearly,  there  is  no 


'ultimate'  meteor  burst  system  but  rather  application 
dependent  design.  With  faster  and  cheaper  hardware  for 


control  and  storage  and  renewed  focus  on  this  complex 
channel  It  is  clear  that  ve  are  In  the  exciting  phase  of 
using  this  natural  phenomenon  for  communication. 
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6  TABLES 


Table  1  -  Order-of-Msgnitude  Estimates  of  the  Pro  partial  of  8pondic  Mataon  (1) 


Mataor 

Maas 

Number  of  This 

Electron  Line 

Radius 

Maas  or  Greater 

Density  (electrons 

Particle* 

(I) 

Swept  Up  by  the 

per  meter  of  trail 

Earth  Each  Day 

length) 

Particlaa  paaa  through  the 
atmosphere  and  faO  to 
the  ground 

D 

8  cm 

10 

— 

Particles  totally  die- 

10* 

4cm 

10* 

integrated  in  the 

10* 

2cm 

10* 

__ 

upper  atmosphere 

10 

0B  cm 

104 

10»* 

1 

0.4  cm 

10* 

10” 

10’* 

0.2  cm 

10* 

10” 

10’f 

0.08  cm 

1°1 

10” 

10° 

0.04  cm 

10* 

10” 

10-4 
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Diurnal  variations  in  meteor  activity. 
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Correlation  patterns  in  the  vicinity  of  a  meteor  burst  receiver 
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Meteor  burt»  geometry.  A  signal  fromtX  is  reflected  lo  RX  by  a  trail 
that  is  h  km  above  the  earth.  TX  and  RX  are  separated  by  a  great  circle 
distance  of  L  km. 


FIGURE  6 


r »’»  *’•  o 


Idealized  time  variations  in  signals  from  (a)  underdense 
and  (b)  overdense  meteor  trails. 
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ABSTRACT 

A  model  l<  described  for  the  enderdenst 
meteor  burst  phenomenon  that  takes  Into  account 
Its  random  decay  time-constant  and  derives  the  CDF 
and  IDF  for  the  burst  duration  as  well  as  the 
Joint  IDF  for  burst  duration  and  decay  time- 
constant.  In  addition  the  time  varying  bit  error 
rate  BER  and  Its  average  Is  calculated  for  binary 
FSK  and  coherent  BISK. 

INTRODUCTION 

The  development  and  understanding  of  a  model 
for  the  behavior  of  a  meteor  burst  Is  crucial  In 
the  study  and  use  of  this  phenomenon  as  a  viable 
coauiunlcatlon  channel,  ref.  (1-7,  11).  For  a 
communication  system  engineer  It  Is  Important  to 
have  an  accurate  understanding  of  the  time  varying 
bit  error  rate.  BER,  characteristics. 

It  Is  known  that  the  SNR  of  the  underdense 
trail  decays  exponential  ly  In  time  with  a  random 
Initial  amplitude  and  a  random  decay  time-constant 
which  In  turn  affect  the  BER  variation  In  time. 
To  date,  however,  Investigators  have  used  either  a 
constant  -  (worst  case)  -  SNR  based  on  fixed  decay 
rate  In  the  calculation  of  the  bit  error  rate, 
ref.  (11,  13)  or  employed  a  SNR  based  on  a 
conditional  expected  value  of  decay  rate  given  a 
fixed  burst  duration,  ref.  (10). 

.This  paper  takes  Into  account  effects  of  the 
random  decay  time-constant  In  Its  derivation  of 
BER  expression.  Using  empirically  reasonable 
distributions  for  Initial  SNR  amplitude  and  decay 
constant  we  derive  the  cumulative  distribution 
function,  CDF,  of  the  burst  duration,  tg,  as  well 
as  the  Joint  PDF  of  decay  time-constant  and  burst 
'duration.  The  SNR  Is  then  expressed  In. terms  of  a 
minimum  SNR  level  n  ,  decay  time-constant, B,  and 
burst  duration,  tg.  This  SNR  expression  Is  used 
In  obtaining  the  appropriate  expression  for  BER 
and  average  BER  In  time  for  th%  two  cases  of 
binary  PSK  and  coherent  BPSK. 

THE  MODEL  FOR  METEOR  BURST 

Meteor  bursts  are  characterised  as  overdense 
or  underdense.  Since  the  underdense  trail 
phenomenon  Is  the  dominant  contributor  to  the 
overall  throughput  we  shall  consider  Its  model, 
ref.  (S). 

•Acknowledgement.  This  work  Is  supported  by 
AFOSR  Grant  No.  77747  and  by  SCS  Telecom,  Inc. 


According  to  ref.  (3)  we  can  write  the  SNR  for 

wndardeme  trail  at: 

(4x)3ni*(l-ttn%cot*a)Nf 

Where, 

P|  Is  the  transmitter  power 

Gg  Is  the  receiver  antenna  gain 

Gf  is  the  transmitter  antenna  gain 
X  Is  the  wavelength  of  the  carrier 
D  Is  the  diffusion  coefficient  of  the 
atawsphere  (7.5  mz/s) 

4  Is  the  angle  of  Incidence  of  the 
transmitted  plane  wave 

0  Is  the  angle  between  the  trail  and  the 
great  circle  path  from  receiver  to 
transmitter 

Rj  Is  the  nominal  distance  between  the 
trail  and  the  stations 

rD  Is  the  nominal  Initial  radius  Of  a 
trail  (0.65m) 

o.  Is  the  effective-echoing  area  of  the 
electron  (ID"28  mZ) 

NF  Is  the  background  noise  floor 
B  Is  the  decay  time-constant  of  the  burst 

a  Is  the  angle  between  Incident  E-vector 

at  trail  and  the  direction  from  the 
trail  to  the  receiver 
0  Is  the  electron  line  density 

Further  simplifications  are  according  to  ref  (10). 
NF  Is  computed  using 

NF  -  kTeW  F, 

when  k  Is  Boltfman's  constant,  T_  Is  (770’k),  H  1$ 
the  receiver  bandwidth  and  Ft  Is  the  external 
noise  figure  -  assumed  to  be  due  to  galactic  noise 
only. 

a  and  P  are  assumed  to  be  unlformally  distributed 
over  (-V/Z,  «/Z)  and  hence  <a  sin  >•<  sinp>  ■  7A . 

Also  It  Is  assumed  ford  (the  angle  of  Incidence 
for  the  transmitted  wave)  that  all  trails  are 
located  at  100km  altitude  and  halfway  between  the 
transmitter  and  receiver. 

For  mathematical  trectabtllty  we  rewrite  eq  (l)  as 
follows: 
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5L  •"<  Ru  tbe  ■  inlawm  aulan  electron  l  ine 
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•»*bet  defines  the  burst  end  (duration)  and 
corresponds  to  t  . 
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<*) 


•q.  ())  yields  an  expression  for  the  Initial  SNA 
level,  SNR(fo), 

«,/» 

1  (4) 


Aq1  • 


He 


which  when  substituted  Into  eq  (2)  provides  us 
with  the  slfnal'to-nolse  ratio  in  terms  of  the 
minimum  SNR  level,  ,  burst  duration,  tB,  and 
decay  ttae  constant,  I. 

(t.  -  t)/B 

SNR(t)  •  n  e  ^  (5) 

CPf  AND  PDF  Of  tt 

Throughout  this  paper  we  shal  1  use  the  PDF  of  q 
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and  by  Integrating  over  I  we  get 
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Me  are  left  now  to  substitute  eq.  (B)  In  eqs. 
(10)  and  (11)  resulting  in  the  Joint  PDF  of  t. 
and  B  (B  -  exponents  I)  ° 
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In  addition  we  shall  assume  the  PDF  of  B  as 
exponentially  distributed 

v,ia)  •  (J-)  t***1  eL  •  .Js  sec  a  »  o 


or  Rayleigh 
f|(») 
where 


•  ,dW 
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I  >  0 
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and  B.  •  .35  sec.  Bi  Is  the  expected  value  of  B 
and  Its  value  used  is  from  the  most  recent 
observations  by  RADC.  (The  most  extensive 
program  for  the  study  of  meteor  burst  phenomenon 
Is  currently  conducted  by  the  Air  Force  (RAOC). 
ref  (6)). 

Using  eq.  (2)  to  represent  SNR  behsvior  In  tim^ 
and  eqs.  (3)  and  (7),  setting  q_  •  q^  't  can  be' 
shown  that  the  Joint  PDF  of  t«  and  B  is  given  by 
«R.  (10) 


and  K | (.)  Is  the  modified  llankel  function  or 
order  1. 

and  the  CDF  of  tB 

Ft|j  (ta)  -  1  -  xK,(x)  tg>0  (,«) 

eqs.  (13)  and  (11)  are  compared  graphically  to 
an  exponential  distribution  with  E(tg)  *  .377 
sec.  and  depicted  In  fig.  1. 

Similarly  we  repeat  the  process  by  substituting 
eq.  (9)  (B  Rayleigh)  Into  eqs-  (10a)  and  (11)  to 
provide  us  with  the  joint  PDF  of  tg  and  B 
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and  CDF  for  t. 
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eqs.  (16)  and  (17)  are  compared  graphically  with 
an  exponential  behavior  where  C(tB)  •  .58  sec. 
and  are  depicted  in  Fig  2. 
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CALCULATION  or  |£R  |t)  Agn  m 


In  general  the  hit  errar  rate  Is  aqual  to  a 
function  g(.)  of  the  SNR: 


BEN  -  f  (SNR  ’  V/R) 


where,  M  Is  the  bandwidth  of  the  receiver  and  R 
Is  the  data  rate,  W/R  •  1.25  as  used  In  the 
CONET  system.  SNR  Is  defined  bjr  equation  (1). 


In  our  discussion  two  cases  are  considered 
binary  FSK  and  coherent  BPS*. 


For  FSK  we  have 
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1  ?5n  •I'.-Urt 
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and  for  BPSK 
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•  2  *rfc  (VT.I5  •  n  e»P  Ut|*t|/(*"b))  (20) 


Multiplying  eqs.  (19)  and  (20)  by  eqs.  (12)  and 
(15)  and  Integrating  over  t.  and  B  yield  the 
average  bit  error  rate  for  FSK  and  BPSK. 


For  FSK:  B  assuwed  Is  exponential  (B^  *  .35) 
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Eqs.  (2la)  and  (21b)  for  •  2  are  deplctad  In 
Fig.  3 


Eqs.  (22a)  and  (22b)  for  •  2  ara  deplctad  In 
Fig.  4 


RESULTS 


Eqs.  (13)  and  (16)  for  the  probabl  lty  dens  I  ties 


der  the  assumptions  that  B  Is 


exponential  and  Rayleigh  respectively  provide  us 
with  the  quantitative  Information  to  evaluate 
events  of  burst  duration  accurately.  Perhaps 
more  Important  their  depiction  suggests  a  good 
aggreement  with  an  exponential  behaviour  for  tg 
-  the  currently  prevailing  dogma  -  particularly 
for  the  case  when  B  Is  Rayleigh.  The  case  for 
which  B  Is  assumed  Rayleigh  seems  correct  for  It 
yields  zero  probability  for  small  B  -  fast 
decays  -  which  aggrees  with  our  physical 
Intuition. 


The  results  for  bit  error  rate  and  average  bit 
error  rate  are  graphed  In  Figs.  3  and  4  and  the. 
expressions  In  eqs.  (21)  and  (22)  can  be 
modified  to  test  the  effects  of  other  modulating 
schemes  and  effects  of  different  average  decay 
time  constants  and  minimum  signal  to  noise 
ratios. 


It  Is  Important  to  note,  however,  that  assuming 
the  same  fixed  maximum  allowable  bit  error  rate 
for  both  modulating  schemes,  FSK  snd  PSK, 
Implies  a  higher  minimum  SNR,  n.  for  FSK  then 
PSK  and  a  shorter  burst  duration  for  FSK 
resulting  In  a  worse  throughput  performance  for 
FSK. 


In  addition,  the  graph  for  BPSK  (the  modulating 
scheme  used  In  the  COMET  System)  shows  lower 
average  bit  error  than  FSK  as  expected.  The 
lower  average  BER  encountered  when  B  assumed 
Rayleigh  as  opposed  to  B  assumed  exponential  can 
be  explained  by  noting  the  Bz  In  the  exponent  of- 
eqs.  21b  and  22b. 


CONCLUSIONS 


A  model  has  been  described  for  a  meteor  burst 
coxmiun  lest  Ion  channel  that  takes  Into  effect  the 
Statistical  variation  In  the  decay  rate  of  the 
slgnal-to-nolse  ratio  of  the  channel.  Both  the 
PDF  for  tg,  the  burst  duration,  and  Joint  POF 
of  tg  and  B,  decay  constant,  were  derived  and 
are  consistent  with  empirical  observations  made 
by  RADC  (6).  The  bit  error  rate  for  two 
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Modulating  schemes  FSK  and  BPSK  were  calculated 
by  averaging  over  the  effect  of  decay  rate  and 
;11ttle  difference  Is  found  whether  the  Initial 
assumption  for  the  distribution  of  the  decay 
tlMe  constant  was  Rayleigh  or  exponential.  The 
expressions  found  for  average  IER  can  be 
Modified  to  test  for  other  coding  techniques. 
The  high  BER  suggests  the  need  for  further  study 
of  the  channel  and  the  use  of  coding  for 
Improving  the  BER. 
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ABSTRACT 


In  this  paper  we  present  the  results  or  an 
analysis  or  a  aetsor-burst  channel  showing  the 
•spec ted  throughput  and  the  probability  or 
ooaplstlng  a  given  aessage  under  two  dirrerent 
scenarios .  In  the  rirst  soenarlo,  the 
transalttsr  Is  presuasd  to  have  knowledge 
(obtained  by  appropriate  probe  algnals)  or  the 
sntlrs  duration  or  a  given  aetsor-burst  (MB). 

In  the  second,  only  the  starting  tlae  or  each  MB 
la  known.  Nuaerlcal  results  baaed  on  the 
analysis  are  provided,  and  it  is  concluded  that, 
through  the  use  of  appropriate  signal  design  and 
error-correction  coding,  sumciently  reliable 
coanunlcations  over  the  MB  channel  is 
achievable. 

I.  INTRODUCTION 
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The  possibility  or  using  ionization  trails 
created  by  aeteors  entering  the  earth's 
ataospher*  to  provide  beyond  llne-or-alght 
coaaunl cat ions  has  been  known  ror  a  long  tlae. 
Several  experlaental  [1],  (2]  and  practical  (33 
aysteas  have  been  built  using  this  aode  of 
propagation.  The  physics  of  the  phenoaena 
generating  these  Ionization  trails  and  the 
resultant  reriection  or  radio  waves  has  been 
described  in  detail  elsewhere  in  the  literature 
til,  t53. 

In  this  paper,  we  examine  the  performance 
of  two  protocols.  A  fixed  length  packet  is 
assuaed.  The  first  protocol  atteapts  to 
optiaally  use  all  available  bursts.  The 
terainsl  to  receive  data  Is  assuaed  to  be 
continually  broadcasting  so  as  to  probe  for  a 
channel  opening.  The  teralnal  sending  the 
message  data  begins  transalttlng  as  toon  as  It 
hears  the  probing  signal.  As  a  result,  the  tlae 
delay  between  the  opening  of  the  channel  and  the 
start  of  the  data  transmission  Is  at  aoat  equal 
to  the  one-way  propagation  tlae.  Once  the 
channel  oloses,  the  probing  signal  disappears 
and  data  transalssion  ceases;  the  search  for 
another  channel  opening  now  begins. 
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The  second  protocol  exaalned  utilizes  the 
trail  duration  less  efficiently.  However,  it  is 
stapler  to  lapleoent  and  It  reduces  the  exposure 
of  terminals  to  detection  by  decreasing  the 
transalssion  requirements.  In  this  protocol, 
the  aessage  sending  teralnal  probes  for  a 
channel.  Channel  openings  are  detected  when  the 
sending  teralnal  receives  a  response  to  Its 
probe.  This  results  in  a  delay  almost  equal  to 
the  two-way  propagation  time  prior  to  the  start 
of  data  transalssion.  However,  for  typical 
trail  durations  (.2  to  1  second),  the  two-way 
propagation  tlae,  which  is  usually  less  than  IS 
nsec,  can  be  considered  negligible  In  aoat 
circumstances.  Using  this  protocol,  only  one 
packet  is  transmitted  per  burst. 

We  assume  packetlzed  data  transmissions 
using  binary  FSK .  Noncoherent  Batched  filter 
detection  Is  used,  and  a  rate  -1/2  Reed-Solomon 
code  Is  employed.  A  similar  analysis  of  the 
expected  performance  or  a  communication  system 
attempting  to  utilize  meteor  bursts  for  packet 
communication  was  presented  in  [5].  However, 
most  of  our  models  dirfer  from  that  of  [5),  and 
we  are  able  to  provide  closed-form  analytical 
solutions  for  three  of  our  four  key 
expressions.  These  are  the  throughput 
expressions  for  each  of  two  modes  of  operation 
of  the  meteor-burst  channel,  as  well  as  the 
probability  of  successfully  completing  a  packet 
for  one  of  the  two  modes.  The  probability  of 
successful  completion  for  the  other  aode  of 
operation  Is  presented  as  an  tnrinlte  sutnatlon, 
similar  to  an  expression  given  In  153. 

II.  ANALYSIS 

In  this  section,  expressions  for  both 
throughput  and  probability  of  correctly 
transmitting  a  given  message  are  presented.  Th* 
supporting  analysis  leading  up  to  these  results 
can  be  found  in  (63.  As  Indicated  above,  two 
different  nodes  of  system  operation  are 
considered.  In  the  first  node,  it  Is  assumed 
that  the  receiver  continuously  transmits  a 
sounder  to  the  transmitter.  Any  time  th* 
transmitter  receives  the  sounder  signal,  It 
knows  that  the  meteor  burst  (MB)  channel  is 
present  and  hence  starts  to  transmit  Its 
message.  As  soon  as  reception  of  the  sounder 
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ll|Ml  onmi,  the  trmulttir  ttopi 
transmitting.  Hi*  trnulMloA  of  •  g  Iran  piektt 
night.  Howe vsr,  span  aany  MB’s,  and  it  la 

that  a  oar  tain  — ount  of  overhead  la 
Hr;  for  tha  transalttar  to  aoqulre  tha  HB 
channel  on  aaeh  and  avary  separate  burst. 
Bpooirioally ,  It  la  assumed  that  t.  seconds  ara 
necessary  for  aaoh  acquisition.  in  aaaanoa 
than,  tha  first  nods  of  oparatlon  oorrsaponda  to 
tha  tranaalttar  knowing  praoiaaly  aaoh  inatant 
of  tlaa  ah  an  tha  MB  ehannal  la  praaant. 

With  thla  sconarlo,  tha  following 
■athenatical  nodal  la  uaad't  Lot  n  bo  tha  total 
nuabsr  of  aataor  bursta  that  ooour  in  aoaa  tlaa 
intarval  T-  saoonda,  and  lot  N  ba  tha  nuaber 
•f  burata  In  T_  saoonda  that  allow  saaotly  J 
pnofceta -per -buret  to  ba  tranaslttad .  than 


P(M.)  -  l  P(Bj |n)P(n)  ,  (1) 

n-l 

uh ara  n  la  taken  to  ba  a  Poisson  randos  variable 
with  probability  aaas  function  given  by 

Id 

Tt  Tn  " 

•  x(r) 

TI  (2) 

nl 

T.  la  the  average  Interval  between 
^Clearly,  thla  nodel  la  only 


p(n) 

In  (2), 
burata. 

approximate,  slnoa  the  aunaatlon  In  (1)  allows 
an  Infinite  nuaber  of  burata  to  occur  In  a 
finite  Intarval  of  tlaa.  However,  if  the 
average  Interval  between  burata  la  ouch  leas 
than  tha  tine  duration  T-,  the  approximation 
should  yield  reasonable  results. 

If  we  define  P  aa  the  probability  that  a 
given  burst  lasta  for  a  duration  equal  to  that 
of  exactly  J  packets,  then 


P(kj|n)- 
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n-N. 
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To  evaluate  P(N  ),  the  distribution  of  the 
duration  of  a  burst  oust  ba  known.  Towards  that 
and,  aasune  the  probability  density  of  a  burst 
duration  la  given  by 


f(db)  -  if-  e  “  ,  (*0 

where  T  la  tha  avarage  length  or  a  burst.  With 
thla  nodal,  it  la  straightforward  to  show 
that  Mj  la  also  Poisson.  That  la. 
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'  To  aa  great  an  extent  as  Is  possible,  we 
use  the  formulation  and  notation  of  [5]. 
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and  there  T.  la  the  duration  of  a  paokat. 

To  obtain  tha  expected  nuaber  of 
lnforaatlon  bits  that  oan  be  tranealtted  over 
the  HB  channel  In  T.  seconds,  whloh  we  refer  to 
as  the  throughput2 .let  «n  be  defined  as 

®a 
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where  I  la  the  nuaber  of  information  bits  per 
packet.  That  la,  N„  la  the  actual  number  of 


"B, 

transnltted  bits  In  one  realisation  of  tha 
process.  Then 
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Consider  now  the  second  protocol,  whereby 
the  transmitter  nakes  no  attenpt  at  nonltoring 
the  presence  or  the  channel  other  than  tha  start 
of  the  burst.  However,  each  tine  a  burst  la 
detected,  an  entire  packet  la  transnltted. 

Then,  with  1  again  equal  to  the  nuaber  of 
Information  bits  per  packet,  the  number  of  bits 
transmitted  across  the  channel  In  a  time 
duration  of  T  seconds  is  given  by 


B2 

where 


1  N 


P’ 


(9) 


Np  Is  defined  as  the  nuaber  of  bursts 
whose  duration  is  greater  than  or  equal 
to  t0*Tp.  Notice  that  the  parameter  t_  Is 
larger  in  this  second  node  of  operation  than  It 
was  in  the  first  mode. 

With  this  model,  the  expected  number  of 
bits  that  can  be  transmitted  over  the  channel 
can  be  shown  to  be  given  by 

i-'T° 
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exp 
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While  the  throughput  expressions  derived 
above  give  some  perspective  regarding  the 
performance  of  a  MB  channel,  at  times  a  more 
meaningful  performance  measure  Is  the 
probability  of  successfully  completing  a  message 
within  a  specified  period  of  time.  Considering 
first  the  mode  of  operation  wherein  the  receiver 


‘  Note  that  this  definition  ignores  the 
effect  of  thermal  noise.  The  throughput 
expressions  given  here  can  be  modified  as  In  [9] 
to  include  errors  due  to  noise  by  assuming  a 
specific  IRQ  model.  However,  noise  effects  are 
Included  In  this  paper  *hen  the  probability  of 
correct  transmission  Is  considered. 
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continuously  sounds  tits  channel,  sup  post  s  given 
packet  is  repeated  for  ss  many  times  ns  possible 
over  s  given  durst Ion  of  tins,  sty  ssoonds. 
Lot  L.  bs  defined  ss  the  greatest  Integer  loss 
than  dr  equal  to  T../T-,  and  let  P_M  be  ths 
probability  of  oorrecll y  completing  the  message3 
on  at  least  one  of  L.  euoceaslve  transmissions 
of  the  sane  packet. 

•If  pp  la  dsflnsd  as  the  probability  of 
asking  anerror  in  a  spaolflc  packet,  given  that 
it  has  been  received  (l.e.,  given  that  the 
otiannal  was  present  long  enough  for  the  packet 
to  be  ooapleted) ,  then  It  oan  be  ahoun  that 

L|  ,  ,  lTe*fltn 

rc*  *  ll  (,-pp  )(n*1)|fr^n*  tT~^ 

n-1  1-1 

Bt0  4  tTP  *  td 

(1*1 )Tp»ntn 
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Note  that  in  deriving  (11),  it  was  assumed  that 
errors  in  all  packets  occur  with  the  sane 
probability.  This  is  not  true  in  reality,  since 
the  strength  of  a  MB  decays  approximately 
exponentially.  However,  If  we  use  the  smallest 
slgnal-to-nolse  (SNR)  in  any  packet  for 
computing  the  probability  of  error  of  all  the 
bits  in  the  packets,  then  (11)  will  be  a  worst- 
ease  result. 

Finally,  for  the  second  protocol.  It  can  be 
shown  that  _  . 
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The  results  derived  up  to  this  point  are  in 
terms  of  the  parameter  Ppt  the  probability  of 
making  an  error  in  the  packet.  Hence  an 
explicit  expression  for  Pp  is  needed  in  order  to 
evaluate  the  overall  system  performance. 

Towards  this  end,  assume  that  binary  FSK  with 
noncoherent  detection  is  the  modulation  foraat 
and  that  the  overall  packet  is  to  be  broken  up 
into  smaller  groups  of  information  digits  which 
are  to  be  encoded  with  a  forward-error 
correction  (FCC)  code  and  then  transmitted  over 
the  MB  ohannel. 

As  indicated  above,  because  the  SNR  of  a 
given  MB  decays  approximately  exponentially  In 
time,  the  probability  of  error  due  to  thermal 
noise  et  the  receiver  will  vary  from  symbol -to- 
syabol.  Therefore,  to  keep  the  analysis 
tractable,  the  probability  of  error  results  will 
be  computed  using  the  SNR  that  is  present  at  the 
end  of  a  MB  (l.e.,  using  the  smallest  SNR  of  the 
burat),  so  that  the  final  results  will  really  be 

^  For  the  purposes  of  this  section,  a 
message  will  be  taken  to  be  a  single  packet. 


upper  bounds.  Also,  again  because  of  the 
exponential  decay  in  SNR,  it  is  reasonable  to 
assume  that  the  channel  will  be  quite  bursty  In 
nature.  Hence  the  PEC  scheme  should  be  one  that 
performs  well  in  the  presence  of  bursty  errors, 
end  for  the  purposes  of  this  paper,  a  Reed- 
Soloson  (RS)  code  will  be  used. 

To  be  specirie,  assume  the  packet  consists 
of  I  bits  or  information,  and  these  I  bits  are 
broken  up  into  J  groups  of  mk  bits  per  group. 

An  (N,K)  RS  oode  la  then  used  to  enoode  each 
group,  where  N  -  2*  -  1.  To  reeeive  ths  packet 
correctly,  all  J  codewords  must  be  correctly 
received.  The  probability  of  this  ooourrlng  Is 

1  -  Pp  -  O  -  P,,)*1  .  03 


(?)  P.‘ 


(1  -  Pm) 


E  Is  the  error -correct Ion  capability  of  the  RS 
code  and  P  is  the  error  probability  of  a  RS 
symbol.  It  Is  given  by 

P„  -  1  *  (1  -  Pj"  .  ( 

9  D 


1  ‘  5  SNR 

pb  *  I  •  <» 

la  the  bit  error  rate  Tor  noncoherent  binary 
FSK. 

111.  NUMERICAL  RESULTS 

Figures  1  and  2  contain  curves  showing  the 
probability  or  completing  a  message  for  the 
first  and  second  protocols,  respectively  (l.e., 
equation  (11)  is  plotted  In  Figure  1  and 
equation  (12)  Is  plotted  in  Figure  2).  The 
abscissa  in  both  cases  Is  Tp,  the  total 
observation  time  during  which  the  message  is  to 
be  completed.  The  curves  In  each  figure  are 
parameterized  by  the  slgnal-to-nolse  ratio, 
vhlch  in  turn  Implies  a  value  for  T,,  the 
average  time  between  bursts.  For  tne  results 
presented  In  Figures  1  and  2,  the  following 
pairs  of  SNR  and  Tj  were  used: 


SNR(dB) 

7 

7.5 

B 

10 

12 


T.  (seconds) 
28 
29 
31 
39 
50 


The  message  duration  (packet  size)  was 
taken  as  2*00  Information  bits,  and  those  2*00 
bits  were  divided  Into  3?  codewords.  Each 
codeword  corresponded  to  a  (31,15)RS  code. 
Finally,  TL  and  Tp  were  each  taken  to  be  0.2 
seconds,  and  tQ  was  taken  to  be  .03  seconds  for 
the  first  protocol  and  .06  seconds  for  the 
second  protocol. 

Pecause  the  value  of  Tj  Increases  as  the 
required  SNR  Increases,  at  some  point  requiring 
a  larger  SNR  becomes  self-defeating,  since  the 


MM*r  or  meteor-bursts  in  a  given  Tp  aooond 
Interval  becomes  too  small  to  maintain  a 
specified  degree  of  reliability.  This  la 
olaarly  indloatod  in  oaoh  of  the  two  riguraa, 
whom  It  la  noon  that,  for  tha  parameter  values 
ohoeen,  an  MR  of  I  dB  provides  tha  best 
moult  a. 

If,  than,  an  8  dB  SMB  thraahold  la  uaad, 
and  If,  for  example,  it  la  daalrad  to  ha *a  a 
probability  of  noaaaga  toaplttlon  of  0.9,  from 
-Figure  1  It  la  aaon  that  tha  observation  ttaa 
auat  bo  about  thraa  alnutaa  whan  tha  flrat 
protocol  la  uaad.  Similarly,  If  tha  aacond 
protocol  la  ohooan,  Plgum  2  lndloataa  that  the 
obaarvatlon  time  ohould  be  increased  to  about 
four  and  one -half  nlnutaa  In  order  to  achieve 
tha  aaaa  performance. 

1 9.  CONCLUSIONS 

An  analysis  of  two  different  protocols  for 
usa  over  an  MB  channel  has  been  presented.  It 
was  found  that  with  sufficient  repetition  of  the 
massage  and  with  a  sufficient  aoount  of  error- 
correction  coding,  reliable  communications  could 
be  achieved.  This  Indicates  that,  at  least  In  a 
back-up  node,  the  MB  channel  Is  quite 
appropriate  for  the  transnlsslon  of  relatively 
short  messages.  In  fact.  It  should  be 
emphasised  that  the  results  presented  here  are 
peaaiwlatlc  results,  In  the  sense  that  they  were 


derived  under  the  assuoptione  that  the 

ins tantsnsexis  SNR  on  the  channel  wee  at  all 

times  equal  to  the  lowest  SNR  that  was  deemed 

acoep table. 
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